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ABSTRACT 

The impending era of wide-field radio surveys has the potential to revolutionize our understanding 
of astrophysical transients. Here we evaluate the prospects of a wide range of planned and hypo¬ 
thetical radio surveys using the properties and volumetric rates of known and hypothetical classes 
of extragalactic synchrotron radio transients (e.g., on- and off-axis gamma-ray bursts, supernovae, 
tidal disruption events, compact object mergers). Utilizing these sources and physically motivated 
considerations we assess the allowed phase-space of radio luminosity and peak timescale for extra¬ 
galactic transients. We also include for the hrst time effects such as redshift evolution of the rates, 
K-corrections, and non-Euclidean luminosity distance, which affect the detection rates of the most 
sensitive surveys. The number of detected events is calculated by means of a Monte Carlo method, 
using the various survey properties (depth, cadence, area) and realistic detection criteria that include 
a cut on the minimum variability of the transients during the survey and an assessment of host galaxy 
contamination. For the detected events we also quantify how well each light curve is characterized 
(e.g., the fraction of sources with a measured rise time, decline time, and peak brightness). We find 
that near-term GHz frequency surveys (ASKAP/VAST, Very Large Array Sky Survey) will detect 
few events: < 30 — 50 on- and off-axis long GRBs and off-axis tidal disruption events, and ^10 — 20 
neutron star binary mergers if 1% of the mergers result in a stable millisecond magnetar. Low- 
frequency surveys (e.g., LOFAR) are unlikely to detect any transients, while a hypothetical large-scale 
mm survey may detect ^ 40 on-axis long GRBs. On the other hand, we find that SKA surveys at 
^ 0.1 — I GHz have the potential to uncover thousands of transients, mainly on- and off-axis long 
GRBs, on-axis short GRBs, off-axis TDEs, and neutron star binary mergers with magnetar remnants. 

Subject headings: gamma-ray bursts: general — radio continuum: general — supernovae: general — 
surveys 


1. INTRODUCTION 

Radio astronomy appears poised for a revolution in 
the study of transient phenomena thanks to the advent 
of wide-field interferometers. At meter wavelengths (tens 
to hundreds of MHz) it is now feasible to image the full 
primary beams of dipole antennae, resulting in new ad¬ 
vanced arrays such as the LOw Frequency ARray (LO¬ 
FAR; van Haarlem et al. 2013), the Murchison Widefield 
Array (MWA; Lonsdale et al. 2009) and the Long Wave¬ 
length Array (LWA; Ellingson et al. 2009). At centimeter 
wavelengths (GHz frequencies) a new generation of wide- 
held facilities is also being developed, such as MeerKAT 
(Booth et al. 2009) and the Australian Square Kilometer 
Array Pathhnder (ASKAP; Johnston et al. 2008a), and 
existing facilities are undergoing upgrades, including the 
Apertif project at the Westerbork Synthesis Radio Tele¬ 
scope (WSRT; Oosterloo et al. 2009) and the expansion 
of the Karl G. Jansky Very Large Array (VLA; Perley 
et al. 2011). The VLA upgrade makes it possible to reach 
greater image depth in a shorter time, raising the pos¬ 
sibility of conducting a new survey combining the sensi¬ 
tivity of the FIRST survey (Becker et al. 1995) with the 
large sky footprint of the NVSS survey (Gondon et al. 
1998). These telescopes and new technologies are ulti- 
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mately paving the way for the Square Kilometer Array 
(SKA; Carilli & Rawlings 2004), which will further trans¬ 
form radio astronomy with its unprecedented scale and 
sensitivity. 

Despite the overall maturity of radio astronomy, sur¬ 
prisingly little is known about what astrophysical sources 
might dominate the transient radio sky. At long wave¬ 
lengths, the brightest sources are expected to result from 
coherent emission processes. Examples include ‘giant 
pulses’ from Galactic pulsars (e.g. Jessner et al. 2005) 
and cyclotron maser emission from brown dwarfs and 
possibly exoplanets (Berger 2002; Hallinan et al. 2008); 
brighter extra-galactic analogs of such events are spec¬ 
ulated to accompany rare violent events, such as giant 
magnetar flares (Lyubarsky 2014) and mergers of binary 
neutron stars (Hansen & Lyutikov 2001), and may be 
connected to the fast radio burst (FRB) phenomenon 
(e.g., Thornton et al. 2013; Spitler et al. 2014). At 
centimeter wavelengths, and on longer timescales, syn¬ 
chrotron sources are expected to dominate the transient 
sky. These are generally powered by the shock interac¬ 
tion between fast ejecta from an energetic explosion and 
dense ambient gas, as in the case of supernovae (SNe) 
and gamma-ray bursts (GRBs). Slow-timescale Galactic 
radio transients are expected to be dominated by a differ¬ 
ent population, such as flares from M dwarfs and X-ray 
binaries, with rates that are believed to be comparable 
to those of extragalactic events (Williams et al. 2013). 

In this work, we consider and study extragalactic syn- 
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chrotron transients. These are the best understood 
events based on radio follow-up of discoveries at other 
wavelengths, they are of proven astrophysical interest, 
and they are amenable to detection by imaging surveys, 
especially those operating at centimeter wavelengths. 
Known examples include core-collapse SNe (Kulkarni 
et al. 1998; Weiler et al. 2002; Berger et al. 2003; Soder- 
berg et al. 2010), GRB afterglows (Frail et al. 1997; 
Chandra & Frail 2012), and relativistic jets from tidal 
disruption events (Giannios & Metzger 2011; Zauderer 
et al. 2011; Berger et al. 2012; Zauderer et al. 2013). 

Searches for these and related types of radio transients 
help to motivate the science cases for all of the telescopes 
and surveys mentioned above. With these next genera¬ 
tion projects coming online in the near future, it is es¬ 
sential to consider how many and which kinds of tran¬ 
sient events will be detected, and how survey strategies 
may be optimized to answer particular science questions. 
Most untargeted surveys to date have detected few if any 
events (e.g., Carilli et al. 2003; Croft et al. 2010a; Bell 
et al. 2011), and despite careful vetting, a substantial 
number of detected candidates have turned out to be 
instrumental artifacts (Gal-Yam et al. 2006; Ofek et al. 
2010; Croft et al. 2011; Frail et al. 2012). Thus, most 
of what is known empirically about the transient radio 
sky consists of upper limits. Complicating matters is the 
fact that these upper limits are often difficult to char¬ 
acterize and interpret, because the expected yield of a 
given survey results from a complex interplay between 
the properties of the transient events in question and 
the survey strategy. While this interplay is often con¬ 
sidered in the derivation of upper limits and event rates 
(Bower et al. 2007; Croft et al. 2011; Mooley et al. 2013; 
Williams et al. 2013), it is challenging to harmonize the 
results from different studies for inter-comparison, or to 
extrapolate them to new surveys. 

While the current results from untargeted surveys are 
of limited use for predicting the yields of the next- 
generation projects, all is not lost since some informa¬ 
tion is available from radio follow-up observations and 
from generic theoretical considerations. The volumetric 
rates of some of the expected events have been explored 
in detail (e.g., Wang & Merritt 2004; Guetta & Della 
Valle 2007; Wanderman & Piran 2010; Kim et al. 2013), 
and targeted radio follow-up of transients triggered from 
other wavelengths has yielded significant insight into the 
radio properties of some of the event classes mentioned 
above. In other cases, detailed numerical models allow 
predictions of radio emission that are well-motivated by 
both observations and theory (e.g., Nakar & Piran 2011; 
van Eerten & MacFadyen 2011). 

Despite this available knowledge, the predicted rates 
of radio transient occurrence and detection — taking 
into account the expected luminosities and durations 
— have generally not received detailed theoretical at¬ 
tention. Frail et al. (2012) estimated event rates ig¬ 
noring cosmological effects such as non-Euclidean lumi¬ 
nosity distances, K-corrections, and time dilation, in¬ 
stead assuming the standard isotropic Euclidean scal¬ 
ing iVan_sky(> Fv) oc where iVaii-sky is the in¬ 

stantaneous number of events on the sky brighter than 
a given flux density, F^,. Other works considered spe¬ 
cific event classes with more detailed analyses (Ghirlanda 


et al. 2014; Feng et al. 2014) but do not treat the surveys 
themselves in detail. As we show in this work, the inter¬ 
play between the transient luminosities and timescales, 
the survey cadence, and the total survey duration is crit¬ 
ical to determining the yield of a given survey, and the 
ability to separate transient from steady sources. 

Here we offer the first detailed and comprehensive as¬ 
sessment of the discovery prospects for extragalactic syn¬ 
chrotron radio transients in the upcoming wide-held sur¬ 
vey era. We do so by considering multiple effects: the 
specihc observing strategies of several proposed and hy¬ 
pothetical surveys, the light curves for a wide range 
of known and theoretical transients spanning a broad 
range of energies, collimation, and ambient densities (and 
hence timescales and luminosities), a judicious choice 
of detectability metrics, and realistic volumetric rates 
and their cosmological evolution. Our goal is to provide 
well-motivated estimates for what various surveys will 
actually discover using the best information currently 
available. Future radio surveys also have the potential 
to discover new classes of transients beyond the ones 
we consider here. However, any transient source pow¬ 
ered by synchrotron radiation from a rapidly expanding 
blastwave obeys fundamental physical constraints, such 
as the relation between peak luminosity and timescale, 
and a characteristic evolution from high to low frequen¬ 
cies. Thus, our investigation of a broad range of the 
synchrotron transient parameter space sets meaningful 
bounds on the possible detection rates even of event 
classes that we do not explicitly consider. 

Our primary outputs in this paper are the numbers and 
basic light curve characteristics of detected transients as 
a function of event class and survey. In addition to this 
information it is also critical to consider what can (and 
cannot) be learned and extracted from each event after 
detection, such as the energy scale, collimation, ambient 
density, and nature of the transient. These will depend 
in part on the ability to robustly measure the peak lu¬ 
minosity, characteristic timescale, and redshift from a 
host galaxy. In a follow-up paper (“Paper H”) we will 
investigate this topic, combining the radio light curves 
produced by our simulations from this paper with syn¬ 
chrotron models and realistic follow-up strategies. Our 
analysis will be conducted with an eye towards optimiz¬ 
ing future transient surveys to yield not just detections 
but also astrophysical insight. 

The paper is organized as follows. We begin by consid¬ 
ering the general characteristics of radio emission from 
extragalactic synchrotron transients (§2.1), and estab¬ 
lishing the particular properties of the events that we 
consider in this paper, including luminosities, timescales, 
and volumetric rates (§2.2). In §3 we describe the char¬ 
acteristics of the radio surveys that we simulate. Our 
Monte Carlo method and detection criteria for the sim¬ 
ulations are described in §4, and their results are sum¬ 
marized in §5. We discuss the implications of our results 
and present our conclusions in §6. 

2. EXTRAGALACTIC TRANSIENTS 

In this work we focus on extragalactic synchrotron- 
emitting transients, powered by shock interaction be¬ 
tween fast ejecta and the ambient medium. The produc¬ 
tion of synchrotron emission in such transients is generic, 
and has been observed from various sources (e.g., GRBs, 
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TABLE 1 

Extragalactic Transient Classes 


Transient 

■R(z = 0 ) 
(Gpc-3 yr-i) 

Ek 

(erg) 

n 

(cm-®) 

di 

TO.15 

(d) 

n.3 

(d) 

7-3 (c) 

(d) 

nso 

(d) 

Ref. 

LGRB, 6 lobs = 0.2 

0.3t 

10 ®i 

1 

1 

320 

62 

31 

4 

1 

LGRB, 0obs = 0.4 

H 

10 ®i 

1 

1 

450 

90 

55 

12 

1 

LGRB, Sobs = 0.8 

4t 

10 ®i 

1 

1 

1900 

230 

150 

150 

1 

LGRB, 0obs = 1.57 

12 t 

10 ®i 

1 

1 

1300 

620 

550 

590 

1 

Low Luminosity LGRB (“LLGRB”) 

500t 

1049 

1 

0.8 

200 

43 

90 

no 

9 

SGRB, Sobs = 0.2 

5t 

IqSO 

10 “® 

1 

220 

110 

90 

no 

6 

SGRB, Sobs = 0.4 

15t 

IqSO 

10 “® 

1 

360 

180 

160 

180 

6 

SGRB, Sobs = 0.8 

60t 

IqSO 

10 “® 

1 

730 

480 

410 

480 

6 

SGRB, e^bs = 1-57 

185t 

IqSO 

10 -® 

1 

1900 

2200 

2000 

650 

6 

On-Axis TDE (“Sw J1644-I-57”) 

0 . 01 * 

1052 

0.1 

1 

3700 

920 

1040 

180 

2-4 

Off-Axis TDE, spherical 

1 * 

1052 

0.1 

0.8 

3700 

900 

900 

900 

5 

NSM: prompt BH 

500* 

3 X 10®° 

0.1 

0.2 

4000 

4000 

4000 

4000 

7 

NSM: stable remnant (“NSM-magnetar”) 

5* 

3 X 10®2 

0.1 

1 

2800 

1300 

1300 

1300 

8 

Type Ib/c SNe (“RSN”) 

5000* 

10^7 

1 

0.2 

870 

120 

55 

1.1 

10 


Note. — ^Scaled with redshift according to star formation rate (Cucciati et al. 2012). ^ Scaled with redshift according to volumetric 
density of supermassive black holes (Sijacki et al. 2014). Light curve duration at observer frequency ly — 0.15 GHz. Light curve 
duration at observer frequency u — 1.3 GHz. Light curve duration at observer frequency u — GHz. Light curve duration at observer 
frequency u — 150 GHz. References: (1) van Eerten et al. 2010; (2) Zauderer et al. 2011; (3) Berger et al. 2012; (4) Zauderer et al. 2013; 
(5) Giannios & Metzger 2011; (6) van Eerten & MacFadyen 2011; (7) Nakar Sz Piran 2011; (8) Metzger Sz Bower 2014; (9) Barniol Duran 
et al. 2014; (10) Soderberg et al. 2008 


SNe, TDEs). Here we consider a wide range of transients 
in terms of their energy scale and ambient density that 
include both relativistic and non-relativistic explosions, 
collimated and spherical outflows, with on-axis and off- 
axis orientations. As a result, these transients span a 
broad range of the luminosity-timescale phase-space that 
is commonly used to map radio transients; see Figure 1. 
We note that some of these transients have been exten¬ 
sively observed in the past using radio follow-up of dis¬ 
coveries at other wavelengths (e.g., on-axis GRBs, SNe), 
some include only a few known examples (e.g., TDEs), 
some are robustly predicted but have not been directly 
observed to date (e.g., binary neutron star mergers, off- 
axis GRBs), and some are purely hypothetical (e.g., neu¬ 
tron star mergers giving rise to a stable millisecond mag- 
netar). While other types of extragalactic synchrotron 
radio transients could be hypothesized, the sources con¬ 
sidered here span a sufficiently broad range of properties 
and volumetric rates to encompass future predictions. 


2.1. General Considerations of Timescale and 
Luminosity for Extragalactic Radio Transients 

For the purpose of placing general constraints on the 
light curves of extragalactic synchrotron radio transients 
we begin by considering the ejection of material with 
a kinetic energy of Ek = 10®^ Ek ,51 erg and an initial 
velocity of Vi = fdiC (corresponding to an initial Lorentz 
factor of Fi = (1 — into an ambient medium 

of constant density n = Iuq cm“^. The ejecta transfer 
their energy to the ambient medium at the characteristic 
radius (i?dec) at which point they have swept up a mass 
comparable to ^ l/F^ of their rest mass: 


Rd 


ec 


/ sek 

\4Trnmpc'^T'^fdf ) 

6 X IQiRm 


( 1 ) 


This occurs at the deceleration timescale: 

tdec«i?dec/2cAr? « Ii5d 

( 2 ) 

Prior to tdec the radio brightness will rise, so this 
timescale defines a minimum peak time for radio tran¬ 
sients (although the peak time could be longer; see be¬ 
low). 

If the observing frequency (r'obs = 1 t'GHz GHz) is lo¬ 
cated above both the synchrotron peak frequency (vm) 
and the self-absorption frequency (zza), then the peak 
brightness is achieved at tdec, and is given by (Nakar & 
Piran 2011): 


TT — nnc.rT,T,rP .,^0.83^1.3 ^0.83 o2.3 n-2 ,,-0.65 

hjz,dec ~ 0.05 mjy eg _]^e^ _ 2 Pj G'l, 28 ^GHz ■ 


( 3 ) 

Here = 10^® Dl ^28 cm is the luminosity distance. We 
also make the standard assumption that electrons are 
accelerated at the shock front into a power-law energy 
distribution N(E) oc E~p with p = 2.3 above a minimum 
Lorentz factor of 7m, and that cb = 0.01 eB ,-2 and Ce = 
0.1 Ce,-! are the fractions of post-shock energy in the 
magnetic field and relativistic electrons, respectively. 

Some of the transients we consider in this work (e.g., 
SNe, binary neutron star mergers) produce only non- 
or mildly-relativistic ejecta and are essentially spherical 
with isotropic emission at all times. Other transients 
(e.g. GRBs, jetted TDEs) produce ultra-relativistic 
ejects in tightly collimated jets® with opening angles of 
6 j ^ 1. However, when viewed in an initial off-axis 
direction, the emission from these sources can also be 
approximated as being spherically symmetric once the 
shocked matter decelerates to sub-relativistic velocities 
and spreads laterally into the observer’s line of sight (e.g., 
Zhang & MacFadyen 2009; Wygoda et al. 2011). At this 


^ GRB jets have F; > 100 and 8j ~ 0.1 (e.g., Frail et al. 2001), 
while TDE jets have Fi 10 and 9j ~ 0.1 (e.g., Metzger et al. 
2012; Berger et al. 2012). 
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point the radio emission is no longer strongly beamed, 
and idee and F^^dec can be approximated using /3i « 1 
(Nakar & Piran '2011). 

From Equations 2 and 3 it is clear that the luminosities 
of synchrotron transients scale linearly in proportion to 
their kinetic energy, but that more energetic events also 

evolve on a slower characteristic timescale, tdec oc 
The timescale can be further increased by cosmological 
time dilation, tdec oc (1 + 2^)1 if such events are energetic 
enough to be detectable at substantial redshifts. 

Explosions in high density environments are also more 
luminous (Equation 3) but exhibit shorter timescales 
(tdec oc High densities could thus in principle 

lead to very luminous short-duration transients. How¬ 
ever, as the density increases synchrotron self-absorption 
becomes important and eventually suppresses the bright¬ 
ness and leads to a longer duration because the relevant 
timescale to reach peak brightness is no longer tdec but 
instead the timescale at which Va = i^ohs, i-O., when the 
optical depth is of order unity. The self-absorption fre¬ 
quency at tdec is given by (Nakar & Piran 2011): 


Va{tdec) ~ 0.8GHz 


1.3 


(4) 


Thus, self-absorption may play a significant role in sup¬ 
pressing the brightness and extending the timescale of 
transients, especially at sub-GHz frequencies. On the 
other hand, the steep self-absorbed spectrum oc 
or i/^/^) can lead to a positive K-correction for sources 
that are energetic enough to be detectable at cosmolog¬ 
ical distances (e.g., on-axis GRB afterglows; Frail et al. 
2006). 

Finally, we note that tdec can be very short for rela¬ 
tivistic sources (tdec oc F^ ®'^^) so in principle such sources 
can exhibit short timescales. However, for relativistic 
sources the relevant peak timescale at radio frequencies 
is not determined by tdec because on this timescale the 
synchrotron peak frequency generically obeys Vm ^ t'obs- 
The value of I'm is determined by the Lorentz factor 
of the minimum energy electrons, jm ~ dOce^-iF (Sari 
et al. 1998), such that for large values of F it is lo¬ 
cated above GHz frequencies until the outflow deceler¬ 
ates significantly. Thus, the radio light curve at t'obs will 
reach peak brightness only when Vm = r'obs (as long as 
i^a ^ i^m] otherwise the peak will occur when i>a = t'obs)- 
The value of Vm for relativistic sources in a constant den¬ 
sity medium (for Vm ^a) is given by (Granot & Sari 
2002 ): 


~ 1.0 GHz ggCg _ 2 t 2 ^ i (5) 

where Fkjso is the isotropic-equivalent energy, which 
takes into account the collimation of relativistic sources, 
the scaling to 10 ^^ erg is used since a typical beaming 
correction of ~ 0.01 will give a fiducial Ek ~ 10 ^^ erg, 
and t = 10012 d is the time since explosion. Thus, rela¬ 
tivistic transients with F^ 3> 1 have a timescale for peak 
emission at radio frequencies that is much longer than 

tdec ■ 

To summarize, for non-relativistic transients, or ini¬ 
tially relativistic off-axis transients, the characteristic 
timescale at radio frequencies is generally determined by 
tdec- However, if the density is sufficiently high, such that 
Va{tdec) ^ i^obs, the peak time will be longer, defined by 


time at which Va = i^obs- Similarly, if the source is rel¬ 
ativistic such that r'm(tdec) ^ t^ohs the peak time will 
also be longer than tdec, defined by the time at which 
i^m = t'obs- This means that for synchrotron transients 
the timescale at a given luminosity cannot be made arbi¬ 
trarily short. The boundary defining the allowed phase- 
space of luminosity and timescale based on the above 
considerations is shown in Figure 1. This boundary is 
based on an initially relativistic source decelerating into 
a constant density ambient medium, and with i>a ^ i^m- 
For any other case (e.g., non-relativistic source, Va ^ Vm) 
the timescale at a given luminosity will be even longer. 

Independent of the generic arguments above, syn¬ 
chrotron sources also obey a basic physical constraint 
that restricts how rapidly they can exhibit significant 
brightness variations, namely their brightness tempera¬ 
ture^ is Tb < 10^^ K (the so-called inverse Gompton 
catastrophe limit; Kellermann & Pauliny-Toth 1969). 
The co-moving brightness temperature of an expanding 
source at a time t is given by (e.g., Kulkarni et al. 1998) 


T'b = 


L, r2 


< 10^2 K, 


( 6 ) 


where « F is the Doppler factor in cases when the 
source is expanding towards the observer. A source can¬ 
not vary on a timescale shorter than its light crossing 
time, 5tmin « t(3/2T'^, where the ~ 1/2F^ factor corrects 
for photon arrival delay in the case of relativistic expan¬ 
sion. Equation 6 can thus be recast as a constrained 
relation between the timescale and luminosity: 




I67r2fcr^i/2r3 


1/2 


80 d r-3/2L:,/3>GAz, 


(7) 

where Fi /,30 = L,y/{10^^ erg s ^ Hz 3) and the right- 
hand term takes Fg = 10^^ K. This leads to the same 
conclusion that more luminous sources have longer char¬ 
acteristic timescales. We show the boundary based on 
the condition Tb < 10^^ K for sources with F/S = 1 
in Eigure 1. This boundary is remarkably similar to the 
one defined by the synchrotron model considerations dis¬ 
cussed above. 

Both considerations imply that mJy-level GHz tran¬ 
sients, which are roughly at the detection threshold of 
planned near-term surveys (Table 2), will have charac¬ 
teristic durations of > 100 d at cosmological distances. 
Shorter durations are possible, but only if the luminosi¬ 
ties are correspondingly lower, making the detection vol¬ 
ume much smaller for such sources; i.e., they will only 
be detectable if they have a high volumetric rate. Eor 
a fixed transient luminosity, timescale, and volumetric 
rate, it is straight-forward to compute the instantaneous 
number of events on the sky above a given flux density 
threshold (F^jim): 


char 


^ / F, 

3 V47rFi.jini 


3/2 


TZtd 


ur 


0.9Ll%ntdur,2(8) 


where TZ is in units of Gpe ^ yr ^ and we use = 1 

mJy. In Figure 1 we show contours of TZ such that one 


^ We note that the variability timescales and brightness temper¬ 
atures for various radio sources were recently compiled by Pietka 
et al. (2015). 
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1.3 GHz 



<5tobs (days) 

Fig. 1.— Phase-space of peak luminosity and timescale for transient synchrotron sources at a fiducial frequency of 1.3 GHz. Sources 
in the region of phase-space above the solid black line will violate the brightness temperature constraint Tg < 10^^ K (Equation 7). An 
independent constraint on the allowed phase space is given by the condition that Um-,^a ^ t^obs (Section 2.1; dot-dashed blue line). The 
green dotted lines indicate the minimum volumetric rate of transients required for at least possible one detection across the entire sky 
for a three-year survey with a limiting sensitivity of 1 mjy (§2.1). Transients with lower peak luminosities (and hence potentially shorter 
peak timescales) can only be detected within smaller volumes and therefore require higher intrinsic rates. These minimum rates can be 
compared to actual transient rates as summarized in Table 1. Finally, colored squares denote the peak luminosity and timescale of the 
various transients considered in this paper: on- and off-axis LGRBs (red), on- and off-axis short SGRBs (gray), neutron star binary mergers 
with prompt black hole formation (NSM: tan), neutron star binary mergers with a magnetar remnant (NSM magnetar: blue), on-axis jetted 
TDEs (SwJ1644+57: light green), low luminosity GRB (LLGRB: brown), and Type Ib/c SNe (RSN: black). We note that all of these 
sources indeed obey the basic luminosity-timescale limits discussed in §2.1. 


event could possibly be detected in a three-year transient 
survey: given values for L^, and fdur we find Vchar(^) = 1 
in Equation 8, adopting tdur = 3 yr if the actual tran¬ 
sient duration is shorter than this value. We find that 
even with this generous definition of detectability, tran¬ 
sients varying on a timescale as short as ^ 3 days are only 
detectable if their volumetric rate is > 10^ Gpc“^ yr“^. 
This is much higher than the rate of all known extra- 
galactic transients with the exception of Type Ib/c SNe 
(although the latter are non-relativistic sources and thus 
vary on much longer timescales than indicated by the 
condition Tb ~ 10^^ K). For more reasonable volumetric 
rates of relativistic transients (< 1 Gpc“^ yr“^; Table 1) 
the minimum timescale is ~ 30 days; this timescale in¬ 
creases to ^ 100 days for a detection rate of a few events 
per year. Thus, extragalactic synchrotron radio tran¬ 
sients with timescales of ^ 100 days are unlikely to be 
detected in upcoming surveys. 


2 .2. Models and Volumetric Rates 

In Figure 2 we show the light curves for the various 
transients that we consider in this paper at four fre¬ 
quencies (0.15, 1.3, 3, and 150 GHz) characteristic of 
planned and hypothetical surveys (Table 2 and §3). As 
can be seen in Figure 1, these transients span a wide 
range in the luminosity-timescale phase-space and there¬ 
fore provide a broad view of extragalactic transients. It 
is also clear from Figure 1 that these transients indeed 
obey the boundaries defined by the basic synchrotron 
model consideration discussed above and by the condi¬ 
tion Tb < 10^^ K. Finally, by comparing the actual volu¬ 
metric rates of these transients (Table 1) to the contours 
showing the rates required for our fiducial level of ob¬ 
servability (one event brighter than ~ 1 mJy on the sky 
over three years), it is already clear that few detections 
will be made in GHz surveys with mJy sensitivity. 

In the following sections we discuss each class of tran- 
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sients in detail. Unless otherwise noted we adopt the 
microphysical shock parameters eg = 0.1, es = 0.01, and 
p = 2.3, typical of those found or used for known extra- 
galactic transients. 

2.2.1. Classical Long GRBs 

The afterglows of long GRBs (LGRBs) are powered by 
the deceleration of an ultra-relativistic jet by interaction 
with the circumburst medium (Meszaros & Rees 1997). 
The isotropic equivalent energies of classical LGRBs are 
large, Uiso ~ 10^^ — 10^"^ erg, but with typical jet open¬ 
ing angles of 9j ^0.1 (Bloom et al. 2003) the beaming- 
corrected energies are ~ 10®° — 10®^ erg (e.g., Frail 
et al. 2001; Berger et al. 2003; Cenko et al. 2012). This 
has been confirmed by late-time radio calorimetry (Frail 
et al. 2000; Berger et al. 2004; Frail et al. 2005; Shivvers 
& Berger 2011). LGRB jets interact with either the in¬ 
terstellar medium or the wind of the massive progenitor 
star, with typical values of n 0.1 — 10 cm“®. 

Here we consider synchrotron emission from the for¬ 
ward shock, neglecting emission from the reverse shock 
which travels into the ejecta at early time (e.g.. Sari & 
Piran 1999). This is justified by the fact that the re¬ 
verse shock emission is generally fainter than the forward 
shock emission at the relevant frequencies of < GHz (e.g., 
Kulkarni et al. 1999; Laskar et al. 2013). However, at the 
high frequencies relevant to millimeter surveys, the re¬ 
verse shock emission can be substantially brighter than 
the forward shock for ^ day (Laskar et al. 2013); we 
discuss the implications of this point in §6. 

The observed light curves of LGRB radio afterglows 
also depend on the observer viewing angle with respect 
to the jet axis, 0obs- We thus consider both on- and off- 
axis LGRB afterglow light curves, based on the relativis¬ 
tic numerical hydrodynamical simulations of van Eerten 
et al. (2010)®, including self-absorption. We utilize a 
model with properties typical of observed LGRBs: 9j = 
9ohs = 0.2, Ek = 10®^ erg, n = 1 cm“®, and p = 2.5, as 
well as off-axis (“orphan”) jets with 0obs = 0.4, 0.8, 1.57 
(Figure 2). On-axis events, though more luminous, rep¬ 
resent only a small fraction, /obs = Jg^sin9d9, of the 
total event rate; the integration from 9i = 0obs/2 to 
^2 = ^obs covers the approximate range of solid angle 
covered by observers with viewing angle 9 k, 9ohs- De¬ 
pending on observing frequency, depth, and cadence of 
the survey, detection rates may be dominated by on- or 
off-axis events. 

For the local volumetric rate of on-axis LGRBs we 
adopt a fiducial value of 77.LGRB(^obs 9-^) = 0.3 
Gpc“® yr“®, motivated as follows. Wanderman & Piran 
(2010) estimate 77lgrb 1 Gpc“® yr“^ for events with 
isotropic luminosity of Ljso ^ 10®° erg s“^ [E-^so 10®^ 
erg for a typical LGRB duration of 10 seconds). The 
LGRB luminosity function measured by Swift is approxi¬ 
mately flat per unit log Ljso below the characteristic value 
Liso ~ 10®^ erg s“^ (corresponding to the same isotropic 
energy ^ 10®° erg in our assumed afterglow calculation 
above). Thus the rate of events with Eiso ~ 10®° erg is 
1/3 of the rate of events with Eiso ^ 10®^ ergs, or « 0.3 
Gpc“® yr“^. For off-axis LGRBs we scale the on-axis 

® http://cosmo.nyu.edu/afterglowllbrary 


rate according to 77.LGRB(^*obs) = 77.LGRB(0obs = ^j)/obs- 
Finally, the rate is assumed to evolve with redshift in 
the same manner as the cosmic star formation rate den¬ 
sity, according to the model of Gucciati et al. (2012). 
Although this may not hold strictly (e.g., Yiiksel et al. 
2008), the uncertainty introduced by this assumption is 
small compared to that from other approximations. 

2.2.2. Low-Luminosity Long GRBs 

A class of low-luminosity LGRBs (LLGRBs) with lower 
total beaming-corrected energies of Ek ^ 10"^° erg and 
only mildly-relativistic ejecta (F^ Ri 2) has been iden¬ 
tified in recent years (Kulkarni et al. 1998; Soderberg 
et al. 2004). LLGRBs appear to be roughly isotropic, 
and it has been suggested that they result from relativis¬ 
tic shock break-out following the core-collapse of massive 
stars with extended envelopes (e.g., Matzner & McKee 
1999, Nakar & Sari 2012). We model the radio emis¬ 
sion from LLGRBs using the model of Barniol Duran 
et al. (2014) as fit to the radio afterglow of GRB 980425 
(Kulkarni et al. 1998). 

We use an LLGRB volumetric rate that is 10 times 
higher than the beaming-corrected LGRB rate (e.g., 
Soderberg et al. 2004; Goward 2005; Guetta & Della 
Valle 2007). Adopting a total on-axis LGRB rate of ~ 1 
Gpc“® yr“^ and a classical LGRB beaming fraction of 
~ 50 we thus take the volumetric rate of LLGRBs to be 
77 -llgrb ~ 500 Gpc“® yr“^. We note that given their 
lower energy scale, LLGRBs are only detectable nearby, 
and therefore redshift evolution of the rate is unimpor¬ 
tant. 

2.2.3. Short GRBs 

Short-duration gamma-ray bursts (SGRBs) are also 
accompanied by afterglow emission, but with beaming- 
corrected energies of Ek ^ 10^° — 10®° erg, lower than 
in the case of LGRBs (Berger 2014). The circumburst 
densities are also theoretically expected and observation- 
ally inferred to be lower than in the case of LGRBs, 
with n oo 0.001 — 0.1 cm“® (e.g., Fong & Berger 2013). 
We therefore adopt the low energy afterglow model of 
van Eerten & MacFadyen (2011) with Ek = 10®° erg, 
n = 10“° cm“°, 9j = 0.2, Cg = cb = 0.1, and p = 2.5. 
These parameters are consistent with the observed on- 
axis optical afterglows of SGRBs (Metzger & Berger 
2012). In addition to the on-axis case {9ohs = 0.2) we 
also consider off-axis events with 9ohs = 0.4, 0.8, 1.57. 

We use a local volumetric rate of on-axis SGRBs of 
"^SGRB = 5 Gpc“° yr“° (e.g., Wanderman & Piran 
2014), and scale it appropriately for off-axis angles as 
described in §2.2.1. As in the case of LLGRBs, redshift 
evolution of the rate is unimportant due to the low lu¬ 
minosity of these events. 

2.2.4. NS-NS Mergers with Prompt Black Hole Formation 

The merger of two neutron stars, or of a neutron star 
and a low mass black hole, results in the ejection of a 
small quantity of mass, Mej ^ 10“° — 10 “^Mq, with 
mildly relativistic velocities, f3i = 0.1 —0.3 (e.g., Rosswog 
et al. 2013; Hotokezaka et al. 2013). The interaction of 
the ejecta with the surrounding interstellar medium will 
result in synchrotron radio emission similar to that of a 
supernova remnant (Nakar & Piran 2011; Hotokezaka & 
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150 MHz 1.3 GHz 



3 GHz 150 GHz 



Fig. 2. — Light curves of the various transients considered in this paper (Table 1) at observing frequencies of 150 MHz (upper left), 1.3 
GHz (upper right), 3 GHz (lower left), and 150 GHz (lower right). Sources considered here include: on- and off-axis long GRB afterglows 
for various observing angles 6 oi,s (LGRB: red); low-luminosity LGRB afterglows (LLGRB; brown); Type Ib/c SNe (RSN; black); on- and 
off-axis short GRB afterglows (SGRB: charcoal); on-axis jetted TDEs (Sw J1644-I-57: green); off-axis TDEs (maroon); neutron star binary 
mergers with prompt black hole formation (NSM; tan) and a magnetar remnant (NSM magnetar; blue). Flux densities are normalized 
to Dl = 10^® cm (z Ri 0.55). In each panel we also annotate the typical flux density of a host galaxy forming stars at a fiducial rate of 
1 M 0 yr“^ (horizontal black bars). See Section 2 for further details. 


Piran 2015). We calculate the radio emission using the 
spherical model of Nakar & Piran (2011) for Pi = 0.2, 
Ek = 3 X 10®° erg (Mej = IO'^Mq), and n = 0.1 cm“®. 

We use a local neutron star merger (NSM) volumetric 
rate of TT-nsm = 500 Gpc“® yr“®, consistent with that 
inferred from the Galactic binary pulsar population (e.g., 
Kim et al. 2013). 

2.2.5. NS-NS Mergers Leaving a Stable NS Remnant & 
White Dwarf Accretion-Induced Collapse 

In some cases neutron star binary mergers may re¬ 
sult in the formation of a long-lived neutron star in¬ 
stead of immediate collapse to a black hole (e.g., Met¬ 
zger et al. 2008; Siegel et al. 2014). This stable rem¬ 
nant is formed rapidly rotating, with a period of P ~ 1 
ms, and correspondingly large rotational kinetic energy 
of Erot ~ 47r^//P^ « 3 X 10®^ erg, where I ~ 10^® g cm^ 
is the NS moment of inertia. If the merger remnant also 
possesses a moderately strong dipole magnetic field of 
B > 10^® G, then its electromagnetic dipole spin-down 
will transfer the rotational energy to the small quan¬ 
tity of ejecta, accelerating it to trans-relativistic speeds. 


Pi ~ Erot/MejC^ > 1, before the ejecta have been de¬ 
celerated by the ISM (Metzger & Bower 2014). A mil¬ 
lisecond neutron star with similar properties may also 
be formed by the accretion-induced collapse (AIG) of a 
rotating white dwarf (Usov 1992). 

We model the radio emission from both neutron star 
mergers leaving a stable remnant and AIC events using 
the Nakar & Piran (2011) spherical blast wave model® 
with Pi = 1, Ek = 3 X 10®^ erg, and n = 0.1 cm“®. The 
energy and velocity of the ejecta are uncertain because 
they depend, among other things, on how efficiently the 
NS outflow couples its energy to the merger ejecta (e.g., 
Bucciantini et al. 2012) and whether a fraction of the 
NS rotational energy is instead lost to gravitational ra¬ 
diation (e.g., Gorsi & Meszaros 2009). The energy will 
also be smaller if the NS collapses to a black hole before 
the bulk of its rotational energy is extracted; however, 

® Piro & Kulkarni (2013) argue that the nascent pulsar wind 
nebula created by AIC will also produce a radio transient, but we 
estimate that the synchrotron radiation from the ISM interaction 
will dominate the nebular emission unless the ISM density is ex¬ 
ceedingly low, n <C cm“®. 
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the much brighter radio emission from stable NS rem¬ 
nants imply that they will dominate the radio detection 
rate if such objects form at all. A higher ISM density 
is expected compared to the case of prompt BH forma¬ 
tion because mergers leaving stable remnants may occur 
preferentially within their host galaxies due to the small 
natal kicks that are expected to accompany the birth of 
the lowest mass neutron stars formed in electron capture 
supernovae (e.g., Belczynski et al. 2008). 

The rate of mergers leaving stable NS remnants is un¬ 
known, as it depends on the equation of state of ultra- 
high density matter and on the mass distribution of bi¬ 
nary neutron stars. We therefore scale the rate to 1% of 
the total NSM rate, i.e., ^ 5 Gpc“^ as we show 

later on, much higher volumetric rates may already be 
ruled out by past radio surveys due to the extremely 
high luminosities of these events (see in particular the 
FIRST/NVSS constraints from Levinson et al. (2002) in 
Figure 3). The rate of AIC is also uncertain, with pop¬ 
ulation synthesis estimates spanning 7?,aic ~ 10—10^ 
Gpc“^ yr“^ (Yungelson & Livio 1998), and the fraction 
of WDs that are rapidly rotating at the time of collapse is 
even less certain. We scale the volumetric rates with red- 
shift according to the cosmic star formation rate density 
(§2.2.1), assuming that the characteristic delay time be¬ 
tween star formation and these events is relatively short, 
<C Gyr (e.g., Leibler & Berger 2010). We note that if 
the typical delays are in fact much longer, the observed 
rate will be suppressed due to the loss of events at ^ > 1. 
Given that both NS mergers with stable magnetar rem¬ 
nants and rotating AIG events are predicted to produce 
similar radio transients and that the rates of both are 
highly uncertain, for notational simplicity we hereafter 
denote both models as “NSM-magnetar.” 

2 . 2 . 6 . On-Axis Jetted Tidal Disruption Event 
(SwJ1644+57) 

The transient event Swift J164449.3 -I- 573451 (here¬ 
after Sw J1644-I-57) exhibited several properties that led 
to its interpretation as a tidal disruption event with a 
relativistic jet (Bloom et al. 2011; Burrows et al. 2011; 
Levan et al. 2011; Zauderer et al. 2011). The relativis¬ 
tic outflow produced long-lived radio synchrotron emis¬ 
sion, which is still detected at the present (Giannios & 
Metzger 2011; Zauderer et al. 2011; Berger et al. 2012; 
Metzger et al. 2012; Zauderer et al. 2013). A second 
TDE candidate with similar X-ray and radio emission, 
Sw J2058-I-05, has been reported by Genko et al. (2012). 

We model the radio emission from on-axis jetted TDEs 
s-t ^ 1 GHz directly using the radio light curves of 
Sw J1644-I-57 as a template (Zauderer et al. 2011; Berger 
et al. 2012; Zauderer et al. 2013). Since data are not yet 
available at 150 MHz (Gendes et al. 2014), we assume 
that the on-axis emission at this frequency is identical to 
the isotropic, off-axis model (see below). This is justi¬ 
fied by the high self-absorption frequency at early times, 
which will strongly suppress the low frequency emission 
for on-axis sources. 

We estimate the local volumetric rate of on-axis jetted 
TDEs in two ways. First, the fact that one event was 
detected by Swift in At « 10 years of monitoring to a 
redshift z = 0.35 (co-moving volume of H « 11 Gpc“^) 
suggests a rate of 77 .swJ 1644+57 l/(HAr) ^ 0.01 Gpc“^ 


yr“^. Alternatively, we can utilize the theoretically ex¬ 
pected TDE rate of ~ 10“® — 10“"^ galaxy"^ yr“^ (e.g., 
Wang & Merritt 2004; Stone & Metzger 2014), which 
given a local galaxy density of ^ 10“^ Mpc“^ corre¬ 
sponds to a TDE rate of ^ 100 — 1000 Gpc“^ yr“^. As¬ 
suming a beaming correction of ~ 100 for Sw J1644-1-57 
(Zauderer et al. 2011; Berger et al. 2012; Metzger et al. 
2012) and that < 10% of TDEs produce relativistic jets 
(Bower et al. 2013; van Velzen et al. 2013), this corre¬ 
sponds to a volumetric rate of <0.1 — 1 Gpc“^ yr“^ 
for on-axis jetted TDEs, consistent with the empirical 
estimate. In what follows we adopt the conservative on- 
axis rate of 77.Swji644-i-57 ~ 0-01 Gpc“^ yr“^ but note 
that the true local rate may be higher by at least an or¬ 
der of magnitude, especially when considering that the 
conditions to produce a mildly-relativistic outflow with 
Fi > 2, which will produce radio emission, may be less 
stringent than those required to produce bright X-ray 
emission (and hence that SwJ1644-|-57 provides only a 
lower limit on the rate). 

We scale the TDE rate with redshift using the volume 
density of ^ 10® — 10 ^Mq black holes, for which we em¬ 
ploy the model of Sijacki et al. (2014) (their Figure 2). 
This model predicts that the TDE rate decreases by a 
factor of ^ 2 between z = 0 and z ~ 4, before decreasing 
more rapidly at z > 4 (although see Hopkins et al. 2007 
who find a more rapid evolution with redshift). Beyond 
the uncertainties in the evolution of the supermassive 
black hole volume density, the TDE rate per galaxy also 
depends on other factors, such as the mass distribution of 
stars in galactic nuclei, which could also evolve strongly 
with redshift. Our assumed rate evolution thus repre¬ 
sents a best guess given the lack of detailed studies. 

2.2.7. Off-Axis Jetted Tidal Disruption Events 

For Sw J1644-I-57 the initial jet Lorentz factor was es¬ 
timated to be Fi ^ 10 (Metzger et al. 2012; Berger et al. 
2012), corresponding to a beaming fraction /b 0.01. 
Off-axis jetted TDEs may thus dominate on-axis events 
at radio frequencies (Giannios & Metzger 2011; van 
Velzen et al. 2011; Donnarumma et al. 2015). Detailed 
predictions for the off-axis radio light curves of jetted 
TDEs are not currently available, but because the emis¬ 
sion is expected to be relatively isotropic and only mildly- 
relativistic at late times (e.g., Berger et al. 2012), we 
assume a spherical blast wave model. We utilize the 
Nakar & Piran (2011) model with Pi = 1, Ek = 10®^ 
erg, n = 0.1 cm“®, Ce = 0.2, cb = 0.01, and p = 2.3, 
chosen to produce a 1.4 GHz peak flux density in agree¬ 
ment with that of Sw J1644-I-57 (since at late times the 
jet has spread laterally and its emission is approximately 
isotropic). 

We assume the same evolution of the volumetric rate 
with redshift as in the on-axis case, but with a local 
volumetric rate of ~ 1 Gpc“® yr“^ which is fff"^ ~ 100 
times the on-axis rate. 

2.2.8. Type Ib/c Supernovae 

Gore-collapse supernovae (SNe) produce radio emis¬ 
sion as the fastest ejecta interact with the dense gas 
of the progenitor’s stellar wind (e.g., Ghevalier 1982). 
Radio emission is observed in both stripped-envelope 
Type Ib/c SNe (e.g., Berger et al. 2003) and Type H 
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SNe (e.g., Weiler et al. 2002). In fact, the few gen¬ 
uine radio transient discoveries to date have been iden¬ 
tified as Type II SNe: a SN in Markarian 297 (Yin 
& Heeschen 1991), FIRST JI2I550.2-kI30654 (Levinson 
et al. 2002; Gal-Yam et al. 2006), and SN 2008iz (Brun- 
thaler et al. 2009, 2010). We note, however, that only 
FIRST J121550.2-1-130654 was discovered in an untar¬ 
geted search (using a comparison of NVSS and FIRST 
with a baseline of about 5 years), while the other two 
events were discovered in targeted observations of start- 
burst galaxies with elevated core-collapse SN rates. The¬ 
ory and observations indicate that Type II SNe are ac¬ 
tually challenging to discover in the untargeted surveys 
we consider here. In particular, these sources are dom¬ 
inated by synchrotron self-absorption and free-free ab¬ 
sorption such that the most luminous events (generally 
Type Iln SNe, with Ln ^ 10^® — 10^® erg s“^ Hz“^) have 
typical timescales of about a decade at ^ 1 GHz (Cheva¬ 
lier 2006), and are thus indistinguishable from steady 
sources; Type II SNe that vary on timescales of year 
(generally Type IIP SNe) have much lower luminosities 
(Lj/ < 10^^ erg s“^ Hz“^) and are thus detectable in a 
negligible volume. This indicates that Type II SNe may 
appear in the source lists of radio surveys, but will not 
be easily identified as transients within the time-frame 
of the surveys. 

Here we instead focus on Type Ib/c SNe, for which 
the radio emission evolves more rapidly due to the pres¬ 
ence of faster ejecta (/? O.I — 0.3; Berger et al. 2003; 
Chevalier & Fransson 2006). The Type Ib/c SNe de¬ 
tected in the radio exhibit a wide spread in peak lumi¬ 
nosity and timescale {L,j^p ~ 10^® — 10^® erg s“^ Hz“^ 
and tp ^ 10 — 1000 d), following the general pattern 
in which the most luminous events evolve most slowly. 
At the present, the radio luminosity function of Type 
Ib/c SNe has not been quantified so here we use data for 
the well-sampled, and relatively typical, event SN 2008D 
(Soderberg et al. 2008), modeled using the formulation 
of Soderberg et al. (2005). We note that the SN model 
is somewhat different from the models we consider above 
because the ejecta have a large spread in velocity, with 
the energy distributed approximately as Ek{>v) ^ 
(Chevalier 1982). 

To estimate the local volumetric rate we note that the 
fraction of Type Ib/c SNe with detectable radio emission 
approaches unity only for the nearest events (< 15 Mpc). 
About ^ 1/4 of events are detected within the distance 
of SN 2008D (« 27 Mpc). We therefore assume a fiducial 
rate of 7?.ibc = 5 x 10® Gpc“® yr“^, corresponding to one 
quarter of the Type Ib/c SN volumetric rate (Guetta & 
Della Valle 2007). 

3. RADIO SURVEYS 

A variety of time-domain radio surveys covering a 
broad range of frequencies, sensitivities, cadences, and 
angular resolutions are slated to come online over the 
course of the next decade and beyond. This section sum¬ 
marizes the set of surveys that we simulate in this work. 
Some are currently in advanced stages of development 
or are even already taking data, while others are still in 
early stages of planning. We also consider hypothetical 
surveys, some of which depend on judicious temporal bin¬ 
ning of planned high-cadence surveys, motivated by the 
transient timescale arguments we presented in §2.1. The 


key parameters of the various surveys are summarized in 
Table 2. 

3.1. ASKAP VAST 

The Australian Square Kilometer Array Pathfinder 
(ASKAP) is a precursor and technology development 
platform for the SKA, under development in Western 
Australia. Its wide field of view (~ 30 deg^) and moder¬ 
ately high sensitivity (^^ mJy beam“® for a 10 s integra¬ 
tion) will enable fast, sensitive all-sky surveys (Johnston 
et al. 2008b). The ASKAP Survey for Variables and Slow 
Transients (VAST; Murphy et al. 2013) comprises several 
surveys aimed at detecting transients in the 1.1 —1.4 GHz 
frequency range: ‘VAST-Wide’ will survey lO'^ deg^ per 
day to a Icr rms sensitivity of 0.5 mJy, while VAST- 
Deep Multi-Field (hereafter ‘VAST-Deep’), will survey 
10^ deg^ to 0.05 mJy (Icr) every ^ year. We also consider 
the planned Deep Single Field survey (hereafter ‘VAST- 
Deep-SF’) with 30 deg^ surveyed daily to 0.05 mJy (Icr). 

We assume that both VAST-Wide and VAST-Deep^ 
will last for 3 years, while VAST-Deep-SF is assumed to 
last for 1 year. Although a dynamical cadence has been 
proposed for VAST-Deep, we assume a uniform cadence 
since details of the dynamical cadence have not yet been 
specified. 

For the daily-cadence surveys VAST-Wide and VAST- 
Deep-SF we also consider the possibility of co-adding the 
daily exposures to create deeper maps with a slower effec¬ 
tive cadence, motivated by the long durations of extra- 
galactic synchrotron transient sources. We construct two 
separate mock surveys (‘VAST-Wide-Stack’ and ‘VAST- 
Deep-SF-Stack’) by averaging 30 separate daily obser¬ 
vations of VAST-Wide and VAST-Deep-SF to achieve 
effective sensitivities of 0.09 mJy and 0.01 mJy, respec¬ 
tively. We note that the latter is limited to the estimated 
confusion limit of ASKAP. 

3.2. LOFAR Radio Sky Monitor 

The LOw-Frequency ARrray (LOFAR) operates be¬ 
tween about 30 and 240 MHz (van Haarlem et al. 2013). 
The 32 core stations of LOFAR will ultimately be able 
to return up to 24 individual beams, which will cover a 
patch of the sky as a Radio Sky Monitor (RSM; Fender 
et al. 2006; Best & LOFAR-UK Consortium 2008; Fender 
2012). The currently operational LOFAR RSM Zenith 
Monitoring Program covers ss 1500 deg^ at « 150 MHz 
to a Icr rms sensitivity of 12 mJy (J. Broderick, pri¬ 
vate communication). A mean cadence timescale of ^ 1 
month is planned, corresponding to 36 observing epochs 
for an assumed three year survey. Beyond the current 
plan, we also consider a more ambitious hypothetical 
survey (‘LOFAR-Expanded’), which we assume reaches a 
factor of 10 times deeper and covers a significantly larger 
sky area (~ 10^ deg^). Though ambitious in relation to 
LOFAR’s existing capabilities, such a program is similar 
to the original RSM. 

3.3. VLA Sky Survey 

^ The terms ‘Deep’ and ‘Wide’ are somewhat misleading; both 
surveys cover the same sky area, with the real distinction being the 
higher cadence and lower per-epoch sensitivity of VAST-Wide. We 
nevertheless adopt this terminology here to avoid confusion with 
existing literature. 
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TABLE 2 

Adopted Radio Survey Parameters 


Name 

^res 

V 

(t(“) 

Area 



Ref. 


(arcsec) 

(GHz) 

(mJy) 

(deg2) 

(d) 



LOFAR 

80 

0.15 

12 

1500 

30 

36 

1 

LOFAR-Expanded 

80 

0.15 

1.2 

lO^' 

30 

36 

1 

SKA-Low 

11 

0.15 

0.004 

10 ^ 

90 

13 

3 

VAST-Wide 

10 

1.3 

0.5 

TiF 

1 

1096 

2 

VAST-Wide-Stack 

10 

1.3 

0.09 

10 "^ 

30 

36 

2 

VAST-Deep 

10 

1.3 

0.05 

lO^' 

365 

4 

2 

VAST-Deep-SF 

10 

1.3 

0.05 

30 

1 

365 

2 

VAST-Deep-SF-Stack 

10 

1.3 

0.01 

30 

30 

36 

2 

SKA 

0.9 

1.3 

0.009 

10 ^ 

90 

13 

3 

SKA-Expanded 

0.9 

1.3 

0.006 

3 X lO'^ 

30 

36 

3 

VLASS-Wide 

2 

3 

0.1 


365 

4 

- 

VLASS-Deep 

0.6 

3 

0.003 

10 

365 

4 

- 

CMB 

840 

150 

1.0 

TiF 

1 

1096 

4 

CMB-Stack-1 

840 

150 

0.31 

10 ^ 

10 

no 

4 

CMB-Stack-2 

840 

150 

0.1 

10 ^ 

100 

11 

4 


Note. — All surveys are assumed to last 3 years, except for VAST-Deep-SF with a duration 
of 1 year. ^“^RMS Icr sensitivity; Cadence; Number of epochs. References: (1) Fender 
2012, J. Broderick (private communication); (2) Murphy et al. 2013; (3) Carilli & Rawlings 
2004; (4) B. Johnson, G. Jones, private communication. 


A series of large surveys with the upgraded Karl G. 
Jansky Very Large Array (Perley et al. 2011), collectively 
known as the VLA Sky Survey (VLASS), has been re¬ 
cently discussed. Although the parameters of the VLASS 
components are far from settled, we simulate several 
large surveys using the most recent information avail¬ 
able to us. We consider surveys operating in the S-band 
(2 — 4 GHz) using a combination of the high-resolution 
A and B configurations. The ‘VLASS-Wide’ survey cov¬ 
ers 10^ deg^ at an rms sensitivity of 0.1 mJy (Icr) over 
4 epochs with a yearly cadence. The ‘VLASS-Deep’ sur¬ 
vey covers a much smaller area of 10 deg^ with an annual 
cadence, achieving a Icr sensitivity of ^ 3/iJy. Although 
it is planned for each epoch of VLASS-Deep to achieve 
the target sensitivity by revisiting the survey footprint 
in 10-30 passes over 3 months, we consider only the 
hnal deep integrations because the timescale for multiple 
passes is comparable to that of extragalactic synchrotron 
transients. Finally, a VLASS ‘All-Sky’ survey covering 
about 34,000 deg^ with a sensitivity of 0.1 mJy (Itr) is 
also being discussed, but it will be comprised of only 2 
epochs, limiting its utility for detailed time-domain stud¬ 
ies. We do not consider this survey here. 

3.4. Millimeter Survey 

In recent years mm and sub-mm telescopes have been 
developed to study the cosmic microwave background 
(CMB), which are sensitive, wide-held, and exception¬ 
ally stable. We consider whether these instruments could 
be employed for transient discovery. We estimate what 
a hypothetical ground-based mm survey could achieve 
by extrapolating the parameters of the proposed design 
for a polarimeter based on lumped-element kinetic in¬ 
ductance detectors (LEKIDs; Araujo et al. 2014). This 
design specihes an aperture diameter of 50 cm, with an 
estimated ^ 3.2 mJy sensitivity (Icr), and a daily sky 
coverage of ~ 10^ deg^ (G. Jones, B. Johnson, private 
communication). For the purpose of our study we con¬ 
sider a large 2-m diameter instrument, achieving 1 mJy 
sensitivity (Icr) with a daily cadence. As in the case of 


ASKAP, we also consider time-averaged versions of the 
data reaching Icr rms sensitivities of 0.3 and 0.1 mJy 
for cadences of 10 days (‘CMB-Stack-1’) and 100 days 
(‘CMB-Stack-2’), respectively. 

One signihcant drawback of the wide-held GMB in¬ 
struments is that they have poor angular resolution (14 
arcmin for the instrument considered here). This lim¬ 
its both the ability to identify non-varying point sources 
due to confusion (e.g., Blain et al. 2002) and the quality 
of event localization. However, the excellent stability of 
CMB telescopes allows the detection of variable sources 
far below the confusion limit, and since synchrotron tran¬ 
sients evolve from high to low frequencies with time, 
follow-up of mm transients with other facilities at cm 
wavelengths could lead to better localizations. 

3.5. Square Kilometer Array 

The planned Square Kilometer Array (SKA) will be the 
most powerful radio telescope ever built (Carilli & Rawl¬ 
ings 2004). The observatory will be located in both Aus¬ 
tralia and South Africa and will initially be composed of 
arrays of three distinct telescope designs denoted SKAl- 
low, SKAl-mid, and SKAl-survey. We consider two hy¬ 
pothetical surveys for transients that could be conducted 
with the SKAl-survey array, which is planned to have 
a field of view of 18 deg^ and to operate in a frequency 
range of 0.65 —1.67 GHz. The first (‘SKA’) assumes that 
100 hours will be spent every quarter in a survey of 10^ 
deg^, resulting in a la rms sensitivity of 9 /rJy. The sec¬ 
ond (‘SKA-Expanded’) assumes that the SKAl-survey 
array will be fully dedicated to a continuous transient 
survey of 3 x 10^ deg^ with a monthly cadence, resulting 
in a sensitivity of 6 /rJy (Ict). This hypothetical survey 
defines a practical upper limit to what untargeted GHz 
transient surveys may be expected to yield in the com¬ 
ing decades. Finally, we consider a low-frequency survey 
at 150 MHz (‘SKA-Low’) assuming observations of 10^ 
deg^ with 100 hours per epoch on a quarterly cadence, 
resulting in a sensitivity of 4 /iJy (Icr). 


4. SIMULATED TRANSIENT SURVEYS 
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In this section we combine the radio light curves and 
volumetric rates of the transient sources described in §2 
with the survey parameters defined in §3 to assess the 
detection rates of various transients, as well as some of 
their basic properties (e.g., redshift distribution), using 
Monte Carlo simulations. 

4.1. Detection Criteria 

Our basic criteria for the detection and identification 
of a transient source are straightforward: (1) a > lOcr 
detection during at least one observing epoch; and (2) 
at least a factor of 2 change in brightness during the 
course of the survey (detection or 5cr upper limit). The 
first criterion is motivated by the signihcant potential for 
spurious signals given the experiences of untargeted ra¬ 
dio transient surveys to date (Gal-Yam et al. 2006; Ofek 
et al. 2010; Croft et al. 2011; Frail et al. 2012), as well 
as by the large number of beams and epochs comprising 
each survey (e.g., ~ 10® — 10^^ beams per epoch for the 
ASKAP, VLA, and SKA surveys). The second criterion 
effectively defines the distinction between sources that 
can be identified as genuine transients during the survey, 
and those that evolve too slowly to be distinguished from 
steady or mildly variable sources. Indeed, some variable 
radio sources can vary by about a factor of two on slow 
timescales potentially leading to contamination if a ro¬ 
bust variability criterion, such as the one we use here, is 
not utilized (Becker et al. 2010; Ofek et al. 2011). 

We note that some previous works on the detectabil¬ 
ity of radio transients have utilized 5 (t as a fiducial de¬ 
tection threshold (e.g., Carilli et al. 2003; Bower et al. 
2007). Such a low threshold will naturally increase the 
predicted number of transient detections. However, we 
do not consider this a realistic threshold since even in 
the idealized case of purely Gaussian noise in wide-held 
radio maps, this will lead to an unwieldy 10® — 10® 
false detections per epoch with ASKAP, VLA, and SKA; 
non-Gaussian effects (e.g.. Frail et al. 2012) may increase 
the number of false detections signihcantly. Moreover, a 
5(7 threshold reduces the dynamic range for detecting the 
actual appearance and/or disappearance of a transient, 
and hence the ability to separate transients from variable 
and steady sources. Finally, the reduced dynamic range 
and low signal-to-noise ratio will also inhibit the ability 
to infer the physical properties and classihcation of the 
transients (explosion time, energy, ambient density, ve¬ 
locity, collimation). Thus, we caution that even if the 
goal of a survey is to simply count the number of radio 
transients (e.g.. Figure 3) the detection threshold should 
be ~ 10(7. 

In principle a transient that does not reach a 10a 
threshold in any single epoch could be detected at some¬ 
what lower significance in multiple epochs and hence be 
recovered. We essentially account for this effect through 
our investigation of the “stacked” surveys (VAST-Wide- 
Stack, VAST-Deep-SF-STAGK, CMB-Stack-1, GMB- 
Stack-2) in which we carry out the transient detection 
in temporally-binned versions of the nominal daily ca¬ 
dence data. On the other hand, for the various surveys 
with a yearly cadence (e.g., VAST-Deep, VLASS), which 
is comparable to or longer than the typical transient du¬ 
rations (§2.1), such time-averaging is not feasible and 
10(7 remains a robust single-epoch threshold to suppress 
a large number of false detections. 


We emphasize that we do not explicitly model the 
complex processes of imaging interferometric data, cata¬ 
loging sources, and measuring fluxes, which are subject 
to a variety of systematics that have proven challenging 
to cope with (e.g., radio-frequency interference, spurious 
sidelobe sources, correlator bugs). Instead we assume 
that these effects can be approximated as a degradation 
of a survey’s limiting flux density and/or effective area. 
Accounting for these systematic effects is part of our mo¬ 
tivation for choosing lOcr as a robust detection threshold. 
We also assume that each survey achieves its character¬ 
istic sensitivity uniformly on the sky, while in practice 
the sensitivity is spatially variable, or, equivalently, the 
effective area is a function of sensitivity (e.g., Groft et al. 
2013; Williams et al. 2013). 

The transient detection criteria in radio surveys should 
also include the effects of contamination by radio emis¬ 
sion from star-forming host galaxies. This has the po¬ 
tential to be more of a problem for radio surveys as com¬ 
pared to optical ones because of the challenges in sub¬ 
tracting sources with complex substructure in interfero¬ 
metric data with sparse uv coverage. Transient searches 
will likely be performed using flux measurements that 
include a contribution from the host galaxy, which will 
reduce the effective change in brightness of a transient. 

Host galaxy contamination is of greatest concern for 
sources that occur exclusively in star-forming galaxies 
(LGRBs, LLGRBs, Type Ib/c SNe), but it will also affect 
some fraction of transients that occur in all galaxy types 
(NS-NS mergers, SGRBs, TDEs). Galaxy contamination 
is also a greater problem for lower luminosity transients 
(e.g., LLGRBs, Type Ib/c SNe) than for those that eas¬ 
ily outshine their host galaxies (e.g., on-axis LGRBs). 
In Figure 2, we mark the host galaxy flux density corre¬ 
sponding to a star formation rate of SFR = 1 Mq yr“®: 

Fi/,gai « 0.3 mJy i^ghz^l ,27 

where we use the mean normalization measured for lo¬ 
cal galaxies (Garilli & Yun 1999) and adopt a spectral 
index of —0.7 relevant for galactic synchrotron emission 
which dominates in this frequency range® (Gondon 1992). 
From Figure 2 it is clear that galaxy contamination will 
affect a large fraction of transient classes at 150 MHz, but 
that only LLGRBs, Type Ib/c SNe, NSMs, and off-axis 
SGRBs will be affected at 1-150 GHz. 

The host galaxy contamination will be mitigated for 
sources that are sufficiently nearby such that their hosts 
will be spatially resolved by the survey interferometer. 
However, with the exception of the VLASS and the SKA 
1.3 GHz surveys, the other surveys have large beam sizes 
of - 10" (VAST and SKA-Low), - 1' (LOFAR), and 
^ 14' (GMB). For a typical galaxy size of a few kpc this 
means that the hosts will only be signihcantly resolved 
for VAST and SKA-Low within only a small distance 
of < 50 Mpc, while for LOFAR and the GMB surveys 
the hosts will generally be unresolved. Taking into ac¬ 
count the transient luminosities, the survey angular res¬ 
olutions, the expected host galaxy brightness at the var- 

® The 150 MHz flux derived this way may be over-estimated by 
a factor of ~ 2 due to the presence of free-free absorption (e.g., 
Williams & Bower 2010). 
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ious frequencies, and the occurrence rate of transients in 
star-forming galaxies we indicate in Tables 3-8 whether 
host galaxy contamination is expected to be a problem 
for transient detection. 

4.2. Monte Carlo Method 

Our method for calculating the detection rates and 
measured properties of the various transients with a 
given survey is as follows. We create a large number 
(> 10^) of mock transient light curves using the models 
in §2 (Figure 2), with random time phasing relative to the 
start of the survey and with random distances, weighted 
by the comoving cosmological volume and, when appro¬ 
priate, star formation rate or SMBH density as described 
in §2. The transient light curves are calculated at each 
redshift, taking into account K-corrections and time di¬ 
lation effects when necessary. The detection significance 
at each epoch is determined for each event and survey 
combination, and the detection criteria are applied to 
determine the sample of detected events. Absolute detec¬ 
tion rates are finally determined by weighting the mock 
sample by the volumetric event rate and the survey du¬ 
ration. A new realization of events is considered for each 
survey separately. For each survey we also quantify the 
fraction of events (/y) that pass the flux threshold crite¬ 
rion but are miseed by the factor of two variability cut. 
This fraction provides a measure of how many transients 
could in principle be recovered with a much longer survey 
duration. 

We note that the issue of detection alone is separate 
from the question of what qualities make a given event 
scientifically useful. We defer a detailed discussion of 
this issue to Paper II. However, some insight into this 
question can be gained based on a few properties that 
are readily available from the mock transient sample. In 
particular, for the transient sources that pass the detec¬ 
tion criteria, we determine the number of epochs with 
> 3 ct detections to assess the number of light curve data 
points. We also determine whether the source is detected 
to rise and/or decline by at least a factor of two relative 
to the measured light curve peak. If a source is only de¬ 
tected to rise or to decline during the survey, then the 
only possible measurement is a lower limit on the peak 
flux and on the rise or decline time, respectively. If the 
transient is detected at greater than half of its maximum 
flux at either the first or last epoch of the survey, then the 
rise or decline time cannot be determined, respectively. 
These operational definitions are important because un¬ 
like in the case of radio follow-up of transients discovered 
at other wavelengths (e.g., GRBs, SNe, TDEs), the ini¬ 
tial time of the event is not known a priori. In Figure 5 
we present the cumulative redshift distributions of the 
transients discovered by the SKA surveys. 

5. SIMULATION RESULTS 

The results of our Monte Carlo simulations are sum¬ 
marized in Tables 3-8 and Figures 3-5. The Tables pro¬ 
vide the number of detected transients for each class and 
survey {N), the mean number of epochs at which a de¬ 
tected transient has a flux density of > 3 (t (h), the mean 
redshift of the detected sources (z), the fractions /rise 
and /fall of detected transients with “measured” rise and 
decline times (as defined in §4.2), respectively, and the 
fraction of detected transients with their peak occurring 


during the survey, /peak- For each transient class and 
survey we also estimate whether galaxy contamination 
will be an issue (‘Yes’, ‘No’, or ‘Maybe’). These designa¬ 
tions are based on the transient peak brightness relative 
to the host galaxy radio emission (for a fiducial SFR = 1 
Mq yr“^), as well as on whether all events in a given 
class are expected to occur in star-forming galaxies. Fig¬ 
ure 4 summarizes the total number of detections for each 
transient class by each survey. 

Before addressing the individual results in detail, we 
consider in Figure 3 the number of sources detected 
across the entire sky at a fixed point in time as a function 
of the survey depth (lOcr) for each radio frequency and 
transient class. The number of sources above a given flux 
density is estimated according to: 


A^all—sky(^ -Ff/) 


n{z')uuz') 

Jo \ + z' 


( 10 ) 


where is the co-moving volumetric density, tdur(-z0 

is the transient duration at redshift z (estimated as the 
timescale over which the flux is greater than one half of 
its peak value), the (1 + z)~^ factor accounts for time 
dilation of the event rate, and the integral is performed 
over the co-moving volume V out to the redshift z(Fi,) 
corresponding to the peak flux density F/.p. We stress 
that since this is a snapshot rate no variability cuts are 
applied, and therefore there is no guarantee that all of 
these sources can actually be identified as transients in a 
real survey®. 

As can be seen in Figure 3 all transient classes fol- 
low the relation A^au-sky ocF^ at high flux densities, 
consistent with the expectation for Euclidean geometry 
and a homogeneous source population. However, sig¬ 
nificant deviations are seen at low flux densities due to 
a combination of K-corrections, non-Euclidean luminos¬ 
ity distance, and cosmological evolution of the source 
population. These effects have not been systematically 
considered in previous studies. For luminous sources, 
A^an-sky(> FJ) is Seen to flatten below a critical flux 
density, corresponding to the depth at which essentially 
all events are visible to the edge of the Universe; for such 
sources deeper searches are ineffective at increasing the 
detection rate. Whether the slope of the distribution 
flattens or steepens approaching this plateau depends on 
the importance and sign of the K-correction. An up¬ 
ward steepening of AiaU-sky(> F^) occurs for sources 
with a strong negative K-correction, which is due to the 
steep positive spectrum of synchrotron self-absorption 
(F/ oc or oc e.g., on-axis LGRBs at < 3 GHz. 

On the other hand, a flattening prior to the plateau oc¬ 
curs when the spectrum is flat or inverted; e.g., on-axis 
LGRBs at 150 GHz. 

Shown for comparison in Figure 3 are available con¬ 
straints from previous radio surveys, as well as from the 
surveys we simulate here. The detection limits are cal¬ 
culated at the 95% confidence level for a fiducial tran¬ 
sient duration of 100 days (as motivated in §2.1). The 


® For instance, at 1.3 GHz, Naii_aky(> G-) is at least an or¬ 
der of magnitude larger for the NSM-magnetar case than all other 
events, but we do not predict a commensurate number of discov¬ 
eries in actual surveys, because ~ 75% of these events do not vary 
substantially over the survey duration (Table 7). 
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150 MHz 1.3 GHz 




3 GHz 150 GHz 




Fig. 3.— Number of sources across the entire sky (eq. [10]) at any time above flux Fu at observer frequency v = 150 MHz (upper left), 
1.3 GHz (upper right), 3 GHz (lower right), 150 GHz (lower right). Sources shown include long GRBs (0obs = 0.2, red; dobs = 0.4, blue; 
dobs = 0.8, brown; dobs = 1.57, green), neutron star merger leaving black hole (olive), NSM-magnetar (light blue), on-axis SGRB (grey), 
off-axis SGRB (orange), on-axis TDE (Sw J1644-I-57; dark green), off-axis TDE (purple). Shown for comparison are representative limits 
placed by the surveys discussed in this paper (Table 2) as well as those placed by past surveys (Appendix A; references: Jaeger et al. 
2012b,de Vries et al. 2004, Croft et al. 2010b, Levinson et al. 2002, Gal-Yam et al. 2006, Ofek et al. 2010, Ofek et al. 2011, Bower et al. 
2007; see also Frail et al. 2012). 


methodology for calculating the upper limits, and a sum¬ 
mary of the results from past surveys are given in Ap¬ 
pendix A. With the possible exception of NSM-magnetar 
events (§2.2.5), which at our fiducial rates approach the 
constraint at 1.4 GHz set by FIRST/NVSS (Levinson 
et al. 2002), none of the previous surveys reach the re¬ 
quired depth/sky area to detect any of the known classes 
of extragalactic transients, consistent with the lack of 
confirmed transients to date (Frail et al. 2012). 

5.1. LGRBs 

The VAST (Deep, Wide, and Wide-Stack), VLASS, 
CMB, and SKA surveys all detect at least one LGRB. 

The wide/shallow surveys VAST-Deep, VAST-Wide- 
Stack, and VLASS-Wide detect ^1 — 2 on-axis LGRB 
afterglows, but this number increases to ^ 15 — 30 for 
off-axis events. The off-axis events are generally detected 
over a few epochs, with a low mean redshifts, z ~ 0.1. 
In addition, ^ 50 — 80% will have measured rise and de¬ 
cline times. These GHz surveys should thus produce a 
reasonably complete sample of LGRB afterglows which 
are located sufficiently nearby to possess detectable SNe 


(detected archivally since the SN optical emission will 
probably have faded by the time of the first radio detec¬ 
tion). However, we note that the SNe may be identified 
first in the optical, in which case targeted radio follow-up 
of such events presents a more profitable approach than 
an untargeted wide-held radio search. 

The wide/shallow surveys detect a greater number of 
LGRBs than the deeper, pencil-beam surveys VLASS- 
Deep and VAST-Deep-SF-Stack, which lead to N 1 — 4 
off-axis events. However, the latter do better than pre¬ 
dicted from a homogeneous source population because 
the LGRB volumetric rate increases by an order of mag¬ 
nitude at the larger redshifts (z > 1) to which these 
surveys are sensitive. 

The stacked GMB surveys detect V ss 40 events with 
an average redshift of z « 0.7. The long effective cadence 
('^ 100 days) of the GMB-Stack-2 survey compared to 
the typical LGRB afterglow duration at 150 GHz renders 
most events detectable at only 1 epoch, but each detec¬ 
tion by the GMB-Stack-1 survey (~ 10 day cadence) will 
occur over ^ 3 epochs. In the latter case, the explo¬ 
sion time of the GRB will typically be constrained to 
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Fig. 4.— Predicted number of event detections {N in Tables 3—8) by survey and event class. Predicted values less than 0.1 are not shown. 
Because of the varying quality of constraints on the underlying rate and light curve models, uncertainties in N vary from one event class 
to another (§2.2). 


< 10 days, sufficiently early to trigger optical follow-up 
and lower-frequency radio observations to confirm the af¬ 
terglow origin and to better localize the event. LGRBs 
are sufficiently rare and bright at z 0.7, that the de¬ 
tection or non-detection of y-rays from the sky location 
and explosion time window could constrain the existence 
of on-axis afterglows that are unaccompanied by a GRB 
(so-called “dirty fireballs”) given the bias towards detect¬ 
ing on-axis events at high observing frequencies. 

The number of on-axis LGRBs detected by millimeter 
surveys may be underestimated since we have only mod¬ 
eled emission from the forward shock, neglecting contri¬ 
butions from the reverse shock which are important at 


early times. At 150 GHz the reverse shock emission could 
be a factor of > 10 times brighter than the forward shock, 
but lasting for only ^ 1 day (e.g. Laskar et al. 2013). Ap¬ 
proximately 8 on-axis TDEs were detected by the 1-day 
cadence GMB survey to redshift z ~ 0.34 considering 
only the forward shock emission. From Figure 2, the 
number of events in the survey with detectable reverse 
shock emission could thus be a factor ^ 20 times larger, 
and extending to much higher redshifts z >2. Alterna¬ 
tively, LGRBs that are detected by the stacked GMB sur¬ 
veys via their forward shock emission could be searched 
at early survey epochs for reverse shock emission, pro¬ 
viding tighter constraints on the time of the GRB and a 
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Fig. 5.— Cumulative number of detected sources with redshift less than 2 for each of the three SKA surveys considered in this paper 
(SKA: black, SKA-Expanded: red, SKA-Low: blue). Different source classes are shown in different panels and with different line styles as 
shown in the legends. Note that the top right panel shows results for both LLGRBs and SNe Ibc. 
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probe of the magnetization of the eiecta (e.g. Giannios 
et al. 2008). 

The large SKA and SKA-Expanded surveys at 1.3 GHz 
detect ~ 700 and 6000 LGRBs, with about 1/3 of 
the events viewed on-axis. On-axis events represent a 
larger fraction as compared to VAST/VLASS since they 
are detected to a higher average redshift, z ^ 1 — 2, 
where the LGRB volumetric rate is an order of magni¬ 
tude larger than locally, while the off-axis events are still 
mainly nearby {z ~ 0.5; Figure 5). SKA-Expanded can 
detect on-axis LGRBs out to z > 3. The light curves of 
most detected sources are reasonably well-sampled, with 
h ^ 3 — 27 epochs, depending on the survey and view¬ 
ing angle. SKA-Low detects 400 LGRBs, mostly far 
off-axis (0obs ^ 7’'/2) and originating from z < 0.3. 

We note that our predicted off-axis LGRB detection 
rates are approximately a factor of ^ 3 times smaller 
than those predicted by Ghirlanda et al. (2014), who 
find that ASKAP VAST should have a detection rate of 

3 X 10“^ deg“^ yr~^, corresponding to ^ 90 events for 
10^ deg^ and a 3 year survey duration. This discrepancy 
results largely from our difference in detection threshold: 
Ghirlanda et al. (2014) adopt a 5cr detection threshold, 
which leads to a factor of ~ 3 times higher predicted 
rate. 

5.2. LLGRBs 

The SKA surveys are the only ones expected to de¬ 
tect LLGRBs (Table 4 and Figure 4). This is due to 
the low luminosity of these events (~ 10^ — 10^ times 
dimmer than off-axis LGRB afterglows), which cannot 
be overcome by their order of magnitude higher volu¬ 
metric rate. SKA-Low and SKA-Expanded detect ~ 20 
and 60 LLGRBs at < 150 Mpc (Figure 5). How¬ 
ever, host galaxy contamination could be crippling in the 
case of SKA-Low due to its relatively poor angular reso¬ 
lution of 11", coupled with the high galaxy flux density 
at 150 MHz. We further note that at the mean distance 
of 150 Mpc, the broad-lined Type Ib/c SNe, which thus 
far appear to accompany all LLGRBs, will be discovered 
independently much earlier in the optical. On the other 
hand, radio-discovered LLGRBs without bright accom¬ 
panying SNe will be of great interest. 

5.3. SGRBs 

At most 2 — 3 SGRBs are detected by the wide/shallow 
surveys including VAST-Wide-Stack, VAST-Deep, and 
VLASS-Wide, with a high probability of being close to 
on-axis (0obs ^ 0.4). No events are detected at larger off- 
axis angles. Much larger numbers, ^ 100, 90, 600 are 
detected by the SKA, SKA-Low, and SKA-Expanded, 
respectively. These events are detected in a modest num¬ 
ber of epochs (n ~ 3 — 10), providing reasonably well- 
constrained rise and decline times. 

The mean redshift of the SGRBs detected by the SKA 
surveys, z ~ 0.2, is sufficiently high that most events will 
not have associated gravitational wave detections with 
Advanced LIGO/Virgo (see also Metzger & Berger 2012, 
who reach a similar conclusion based directly on the local 
SGRB rate). For the few events within the gravitational 
wave detection range, the (albeit poor) sky localizations 
from Advanced LIGO/Virgo will be sufficient to uniquely 
associate the gravitational wave source with the SGRB 


radio afterglow, even given a ~ year uncertainty in coin¬ 
cidence time. 

5.4. NS-NS Mergers with Prompt BH Formation 

VLASS-Wide-Stack, VAST-Deep, and VLASS-Wide 
detect < 0.1 NS-NS mergers, although this number is 
uncertain by at least an order of magnitude given the un¬ 
certain local merger rate. These numbers are suppressed 
by an order of magnitude because a substantial fraction 
(/v ^ 0.8 —0.9) of nominal detections are lost to the vari¬ 
ability cut, due to the slow evolution of NS-NS merger 
light curves relative to the survey durations (Figure 2). 
SKA and SKA-Expanded detects ~ 2 and ~ 10 events, 
respectively, with z ^ 0.02, while SKA-Low detects ~ 50 
with z < 0.05. All detected events are sufficiently close 
for Advanced LlGO/Virgo to detect their gravitational 
wave emission, which will happen well in advance of the 
radio light curve maxima. 

5.5. NS-NS Mergers with Stable Magnetar Remnants & 
White Dwarf AIG 

NSM-magnetar events are among of the most promis¬ 
ing transient classes due to their high luminosities. The 
wide/shallow surveys VAST-Wide-Stack, VAST-Deep, 
and VLASS-Wide detect ^ 10, 20, 5 events, respectively 
to an average redshift of z ^ 0.2 — 0.4. Due to the long 
timescales of these events, almost all detections occur at 
every epoch, and the peak is measured in only ^ 30—60% 
of the cases. SKA, SKA-Expanded, and SKA-Low detect 
200, 1000, 1200 events, respectively, to much higher 
redshifts of z ^ 0.8 — 1.4. Unlike the shorter-wavelength 
surveys, SKA-Low has the capability to detect NSM- 
magnetar events to z > 2 (Figure 5). We stress that 
while the optimistic detection rates make these sources 
of great interest for future radio surveys, the fraction of 
mergers that lead to stable magnetars (as well as the 
overall NS-NS merger rate) is highly uncertain. Upcom¬ 
ing surveys will constrain the true occurrence rate. 

NSM-magnetar represent the kind of source for which 
untargeted radio transient surveys have more discovery 
potential (Metzger & Bower 2014) because they are ex¬ 
pected to be radio-luminous without having counterparts 
that are more readily discovered at other wavelengths. 
These events may also be accompanied by luminous X- 
ray and optical emission following the merger or AIG 
(e.g., Metzger et al. 2008; Yu et al. 2013; Metzger & Piro 
2014), but the short duration and low luminosities of 
these putative emission processes will make their discov¬ 
ery challenging as compared to at radio frequencies. 

5.6. On-Axis Jetted TDEs (Sw J1644+57) 

CMB-Stack-2 detects ~ 2 on-axis TDE, with z ~ 1, 
while the wide/shallow GHz surveys VAST-Wide-Stack, 
VAST-Deep and VLASS-Wide detect < 1 on-axis TDEs 
(Table 8). While the peak flux density of these events is 
higher at 3 GHz than at 1.3 GHz, VLASS loses a greater 
number of events to the variability cut (/v ~ 0.5) as 
compared to VAST (/y ^ 0.3 — 0.4). The pencil beam 
survey VLASS-Deep detects fewer events than would be 
predicted even assuming a Euclidean geometry and ho¬ 
mogeneous source population because on-axis TDEs are 
so luminous that their detection rate is limited not by 
the flux density threshold of the survey, but by the de¬ 
creasing SMBH density at z > 4 and the greater loss of 
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events due to the variability cut introduced by time di¬ 
lation (/v = 0.75). The same effect limits the number of 
detections by the various SKA surveys to < few, despite 
their much greater sensitivity. In essence, on-axis jetted 
TDEs are rare events, even if one can probe the entire 
observable universe. 

The long durations of jetted TDE afterglows result in 
most events being detected in most survey epochs, e.g., 
h « 7 for CMB-Stack-2, yet the duration is sufficiently 
short that most will have measured rise times. Although 
the initial epoch of jet launching will only be constrained 
by radio observations to within ~ 100 days, the rarity of 
jetted TDEs implies a good chance to unambiguously as¬ 
sociate the radio sources with prompt gamma-ray emis¬ 
sion detected by Swift or Fermi, either via a previous 
trigger or through an archival search (e.g., Cenko et al. 
2012). The soft X-ray emission might also still be de¬ 
tectable, as luminous X-rays accompanied Sw J1644-I-57 
for about 500 days after the TDE (Zauderer et al. 2013). 
We also note that the radio detection rates could be 
higher than our predictions if not all radio-producing on- 
axis TDEs are accompanied by luminous X-ray emission, 
since our rates are derived from X-ray-discovered events. 
However, the typically high redshift of detected events of 
z F 0.6 implies that constraining the transient location 
to the nucleus of the host galaxy with sufficient accu¬ 
racy to confirm a TDE origin will be challenging, likely 
requiring follow-up JVLA or VLBA imaging. 

5.7. Off-Axis Jetted TDE 

Off-axis jetted TDEs exhibit higher detection rates, 
with « 13, 20, 6 events detected by VAST-Wide-Stack, 
VAST-Deep and VLASS-Wide, respectively (Table 8). 
These events occur at lower average redshift (z ~ 0.2 — 
0.3; Figure 5) than the on-axis events and are also de¬ 
tected at most epochs of the survey (h « A^obs)- Detected 
events are reasonably well-characterized, with rise and 
decline times determined for ^ 70% and ^ 30% of events, 
respectively. The typically large viewing angle with re¬ 
spect to the jet axis and relatively small typical distances 
imply that VLBI imaging of off-axis TDEs could allow 
the extended jet structure to be resolved (Giannios & 
Metzger 2011; Mimica et al. 2015), potentially providing 
an independent age for the system. 

LOFAR-Extended detects less than one event [N « 
0.2) due in part to the large number of events lost to the 
variability cut (/v ~ 0.9). The millimeter surveys also 
detect < 1 off-axis TDEs due to the suppressed peak flux 
at high frequencies as compared to the on-axis case. 

We note that our overall TDE detection rates are much 
lower than those predicted by Frail et al. (2012), primar¬ 
ily because these authors assumed a beaming correction 
to the SwJ1644-|-57 rate of 10^ (instead of 10^ here). 
They also assumed that the off-axis TDE light curve at 
5 GHz is identical to the on-axis case (Sw J1644-I-57), 
whereas at 3 GHz we find that the peak flux density is a 
factor of 6 times dimmer in the off-axis case (Figure 2). 
Our assumptions result in a predicted TDE radio detec¬ 
tion rate that is two orders of magnitude lower than in 
Frail et al. (2012). 

5.8. Type Ib/c SNe 

No Type Ib/c supernovae are detected in any of the 
surveys (Table 4), except perhaps the SKA surveys. 


which may detect a few events. As in the case of LL- 
GRBs, these low detection rates result from the low 
brightness of these events (> 10^ times dimmer than 
off-axis LGRB afterglow), which cannot be overcome by 
their higher volumetric rates. Also, as in the case of LL- 
GRBs, host galaxy contamination is a significant prob¬ 
lem, especially for SKA-Low. 

Although it has been argued that perhaps half of 
all core-collapse SNe are dust obscured (Horiuchi et al. 
2011), the actual fraction will need to be much higher 
than this (^ 90%, i.e., a true core-collapse SN rate that is 
10 times higher than inferred from optical surveys) to al¬ 
low a meaningful rate measurement from untargeted ra¬ 
dio surveys due to the small number of detectable events. 
Because the low-luminosity events do not evolve rapidly 
(typical timescales > 100 days at 1.4 GHz) and are only 
detectable in the local universe, targeted searches of lo¬ 
cal galaxies with a few epochs per year are likely to be 
much more efficient for investigating the prevalence of 
obscured event populations (§2.2.8). 

5.9. Summary 

Figure 4 graphically summarizes the expected number 
of event detections for each source class and survey. We 
briefly summarize the results of our simulations as fol¬ 
lows: 

• No transients are expected to be detected with 
LOFAR (even the hypothetical LOFAR-Expanded 
survey) due to a mismatch between the survey sen¬ 
sitivity and the low peak luminosities of transients 
at sub-GHz frequencies (self-absorption). The long 
durations at these frequencies also lead to many 
events appearing as steady sources. 

• At GHz frequencies, the wide surveys VAST-Wide- 

Stack, VAST-Deep, and VLASS-Wide are the only 
pre-SKA surveys producing more than a handful 
of detections, including ~ 15 — 30 classical LGRBs 
(mostly off-axis), 10 — 40 TDEs (mostly off-axis), 

^ 2 on-axis SGRBs, and ^5 — 20 NSM-magnetar 
events. The deep pencil-beam surveys may detect 
a few transients. However, even though they are 
sensitive to events at higher redshifts, where the 
volumetric rate may be higher (e.g., in the case 
of LGRBs) and positive K-corrections become im¬ 
portant, this is insufficient to overcome their lower 
etendue compared to the wide surveys. 

• The GMB-Stack surveys at 150 GHz detect only 
on-axis events: ^ 40 LGRBs and ~ 2 jetted TDEs. 

• The SKA surveys (in particular our hypothetical 
best-case SKA-Expanded survey) detect hundreds 
or thousands of LGRBs (on- and off-axis), off-axis 
jetted TDEs, and NSM-magnetar events. They 
also detect smaller numbers ('^ 10 — 100) of LL- 
GRBs, SGRBs, and NS-NS mergers with prompt 
black hole formation. 

• Essentially none of the surveys (with the excep¬ 
tion of a few detections with SKA-Expanded) de¬ 
tect Type Ib/c SNe due to their low luminosity. 
While core-collapse SNe have been discovered in 
radio observations, these events are not efficiently 
discovered in untargeted surveys. 
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6 . CONCLUSIONS 

We presented general arguments about the con¬ 
strained timescale-luminosity relation of extragalactic 
synchrotron transients, and utilized the light curves for 
a wide range of known and hypothetical transients and 
the properties of multiple planned and hypothetical ra¬ 
dio surveys spanning 150 MHz to 150 GHz to explore the 
detection rates and basic properties of radio transients. 
Our finding can be summarized as follows: 

• The basic physics of synchrotron-emitting extra- 
galactic transients place significant constraints on 
their possible radio light curves. Of particular rel¬ 
evance to radio transient surveys are the facts that 
luminous events cannot evolve on arbitrarily short 
timescales, and that the emission evolves from high 
to low frequencies. Events with radio luminosities 
sufficient to be detected in appreciable numbers 
with future radio surveys (e.g., LGRBs, TDEs) 
have characteristic timescales of > 100 days at GHz 
frequencies and about an order of magnitude longer 
at - 0.1 GHz. 

• Our Monte Carlo simulations predict detection 
rates for planned and hypothetical radio transient 
surveys, for the first time using an end-to-end 
model including realistic volumetric rates and light 
curves, and simulated survey strategies. We ac¬ 
count for all relevant cosmological effects, includ¬ 
ing time dilation, K-corrections, rate evolution, and 
non-Euclidean luminosity distance, which are im¬ 
portant in the context of the most sensitive surveys 
such as SKA. 

• Low-luminosity transients require volumetric rates 
comparable to the local core-collapse SN rate to 
be discoverable in untargeted radio transient sur¬ 
veys. They will therefore account for at most a 
few events in all pre-SKA surveys, and even with 
SKA will number ^ 10 — 100 in the entire survey. 
Such transients are therefore best studied through 
targeted radio follow-up of discoveries from other 
wavelengths, if they produce emission outside of 
the radio band. However, even in the case that they 
produce only radio emission, or are dust-obscured, 
a more profitable approach may be to target nearby 
galaxies. 

• The small number of expected detections in the 
pre-SKA surveys (and the hundreds to thousands 
of detections with SKA) are dominated by rela¬ 
tivistic sources (LGRBs, TDEs, NSM-magnetars) 
both on- and off-axis. These events are detected 
at cosmological distances where redshift determi¬ 
nations from host galaxies may prove challenging, 
and the ability to locate them to specific regions of 
the hosts (e.g., the nucleus in the case of TDEs) 
may be limited. 

• Our hypothetical millimeter surveys will detect 
mainly on-axis LGRBs and TDEs. The detection 
rates we calculate may be enhanced if additional 
events are detected by their reverse shock emis¬ 
sion, which is significantly brighter than the for¬ 
ward shock emission but shorter lived. LGRBs 


discovered by the stacked GMB surveys via their 
forward shock emission could also be searched at 
early survey epochs for reverse shock emission. 

• Our predicted detection rates are generally lower 
than previous claims (Figure 4). We predict that 
the most successful pre-SKA survey will be VAST- 
Deep, with ^ 70 extragalactic transient detections 
across all considered classes. In some cases, the 
lower numbers that we predict are due to more con¬ 
servative assumptions about the underlying events 
(e.g., §5.7). Another factor is our lOcr detection 
criterion, which is generally higher than previous 
studies, but also more realistic given the scale of the 
surveys and lessons about false-positive candidates 
from past searches (Frail et al. 2012). Our simula¬ 
tions also highlight the impact of slow evolutionary 
timescales at low frequencies which will restrict the 
ability to separate transient sources from variable 
and steady sources. We also note that some pre¬ 
vious studies have overestimated the effective ar¬ 
eas of surveys that observe at cadences signihcantly 
higher than the relevant evolution times. 

• We find that host galaxy contamination is not ex¬ 
pected to be a problem for most detected tran¬ 
sients, because the detections are dominated by 
luminous events that outshine their host galax¬ 
ies for moderate star formation rates of ^1- 
10 Mq yr“^. The small numbers of detections of 
the low-luminosity transient classes will be reduced 
at < GHz frequencies, however, if the host galaxy 
emission cannot be subtracted cleanly. This is es¬ 
pecially relevant for SKA-Low, because at low fre¬ 
quency only a few radio transients (on-axis SGRBs, 
NSM-magnetars, TDEs) manage to signihcantly 
outshine their host galaxies, and the resolution of 
low frequency transient surveys is relatively poor 
(11" for SKA-Low). 

• For near-term surveys, we hnd that areal den¬ 
sities of realistic extragalactic radio transients 
generally follow a Euclidean scaling relation, 

A^aii—sky(^A’t/) oc Fy ^ , where A^aii—sky(^A'iy) is 
number of transients on the sky brighter than some 
hux density cutoff, Fy (Figure 3). As a rule of 
thumb, non-Euclidean effects are not signihcant for 
transients brighter than ~1 mJy in almost all cir¬ 
cumstances, and for some sources they can be ig¬ 
nored down to flux densities of ^ 10 /rJy. However, 
for the more sensitive SKA surveys, which are the 
only surveys expected to detect a large number of 
transients, non-Euclidan effects and cosmological 
source evolution cannot be neglected. 

• The small number of predicted events in pre-SKA 
surveys will inhibit attempts to measure event rates 
and beaming fractions, since the uncertainties in¬ 
duced by small-number statistics will be signihcant. 
It is also unclear if the various transient sources can 
be distinguished from each other in surveys that 
span only a small number of epochs. 

• The SKA will detect hundreds to thousands of 
events. These will be dominated by LGRBs but 
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will also include SGRBs (primarily on-axis), off- 
axis TDEs, and NSM-magnetar events. 

• NSM-magnetar events have the potential to be a 
major source of radio transients and may dominate 
the yield of SKA-Low. However, the underlying 
rates and light curves are currently highly uncer¬ 
tain; upcoming surveys will constrain them. 

We have not considered the question of what astro- 
physical information can be extracted once a radio tran¬ 
sient is discovered. As we have emphasized, the gen¬ 
eral lag of radio emission behind that at other bands 
makes the detailed characterization of radio-discovered 
transients challenging. Paper II will investigate this topic 
in detail. 
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support from NASA Fermi grant NNX14AQ68G, NSF 
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TABLE 3 
LGRBs 


Survey 

N°- 



fd 

Jy 

rise 

•/fall 

rg 

•'peak 

gal?'^ 

N 

n 

z 

/v 

/rise 

/fall 

/peak 

gal? 




^obs 

= 0.2 






^obs 

= 0.4 






LOFAR 

0.0 

26(37) 

0.004 

0.48 

0.75 

0.25 

0.58 

M 

0.0 

25(37) 

0.004 

0.48 

0.82 

0.18 

0.61 

M 

LOFAR-E 

0.0 

26(37) 

0.012 

0.46 

0.77 

0.23 

0.61 

M 

0.0 

25(37) 

0.012 

0.49 

0.82 

0.18 

0.60 

M 

SKA-Low 

5.1 

10(13) 

0.249 

0.28 

0.68 

0.54 

0.64 

M 

12 

9(13) 

0.243 

0.38 

0.75 

0.28 

0.66 

M 

VAST-W 

0.09 

217(1096) 

0.074 

0.00 

0.91 

0.90 

0.91 

N 

0.18 

277(1096) 

0.062 

0.00 

0.87 

0.87 

0.88 

N 

VAST-W-S 

1.4 

7(37) 

0.184 

0.00 

0.91 

0.90 

0.88 

N 

2.6 

10(37) 

0.152 

0.00 

0.87 

0.87 

0.85 

N 

VAST-D 

1.4 

1.1(4) 

0.219 

0.01 

0.73 

0.76 

0.00 

N 

3.6 

1.2(4) 

0.184 

0.01 

0.75 

0.74 

0.00 

N 

VAST-D-SF 

0.0 

170(366) 

0.252 

0.00 

0.75 

0.76 

0.78 

N 

0.010 

208(366) 

0.203 

0.00 

0.68 

0 .68) 

0.71 

N 

VAST-D-SF-S 

0.16 

6(13) 

0.847 

0.00 

0.75 

0.73 

0.73 

N 

0.16 

8(13) 

0.545 

0.00 

0.67 

0.66 

0.65 

N 

SKA 

140 

3(13) 

0.914 

0.00 

0.89 

0.90 

0.31 

N 

150 

3(13) 

0.589 

0.00 

0.85 

0.84 

0.64 

N 

SKA-E 

2060 

8(37) 

1.822 

0.00 

0.89 

0.87 

0.87 

N 

1300 

10(37) 

0.826 

0.00 

0.85 

0.84 

0.84 

N 

VLASS-W 

1.9 

1.0(4) 

0.310 

0.00 

0.76 

0.75 

0.00 

N 

4.2 

1.0(4) 

0.223 

0.00 

0.75 

0.75 

0.00 

N 

VLASS-D 

0.45 

1.0(4) 

1.892 

0.00 

0.76 

0.75 

0.00 

N 

2.0 

1.1(4) 

2.066 

0.00 

0.75 

0.74 

0.00 

N 

GMB 

8.2 

25(1096) 

0.343 

0.00 

0.99 

0.99 

0.99 

N 

1.1 

61(1096) 

0.117 

0.00 

0.98 

0.97 

0.98 

N 

CMB-Sl 

36 

2.9(110) 

0.568 

0.00 

0.99 

0.99 

0.05 

N 

5.6 

7(110) 

0.199 

0.00 

0.97 

0.97 

0.58 

N 

GMB-S2 

32 

1 .0(11) 

0.847 

0.00 

0.90 

0.91 

0.00 

N 

8.4 

1 (11) 

0.282 

0.00 

0.90 

0.91 

0.00 

N 




^obs 

= 0.8 






^obs 

= 1.57 






LOFAR 

0.0 

29(37) 

0.004 

0.44 

0.60 

0.40 

0.34 

M 

0.0 

32(37) 

0.004 

0.18 

0.50 

0.50 

0.37 

M 

LOFAR-E 

0.008 

29(37) 

0.011 

0.44 

0.61 

0.39 

0.36 

M 

0.04 

32(37) 

0.014 

0.19 

0.46 

0.54 

0.34 

M 

SKA-Low 

45 

10(13) 

0.203 

0.49 

0.72 

0.28 

0.43 

M 

320 

12(13) 

0.252 

0.23 

0.47 

0.53 

0.28 

M 

VAST-W 

0.40 

439(1096) 

0.050 

0.00 

0.79 

0.82 

0.86 

N 

0.42 

743(1096) 

0.032 

0.00 

0.62 

0.64 

0.75 

N 

VAST-W-S 

6.2 

15(37) 

0.121 

0.00 

0.78 

0.81 

0.78 

N 

5.4 

26(37) 

0.073 

0.00 

0.57 

0.63 

0.60 

N 

VAST-D 

11.8 

1.5(4) 

0.149 

0.00 

0.73 

0.74 

0.01 

N 

13.0 

2.5(4) 

0.097 

0.00 

0.62 

0.57 

0.25 

N 

VAST-D-SF 

0.50 

435(1096) 

0.166 

0.00 

0.79 

0.81 

0.85 

N 

0.04 

769(1096) 

0.099 

0.00 

0.60 

0.65 

0.75 

N 

VAST-D-SF-S 

0.68 

15(37) 

0.404 

0.00 

0.79 

0.78 

0.78 

N 

0.42 

26(37) 

0.205 

0.00 

0.58 

0.58) 

0.57 

N 

SKA 

270 

5(13) 

0.421 

0.00 

0.77 

0.78 

0.71 

N 

160 

9(13) 

0.209 

0.00 

0.58 

0.59 

0.52 

N 

SKA-E 

1900 

16(37) 

0.567 

0.00 

0.78 

0.78 

0.79 

N 

860 

27(37) 

0.254 

0.00 

0.55 

0.60 

0.54 

N 

VLASS-W 

7.6 

1.1(4) 

0.141 

0.00 

0.75 

0.76 

0.00 

N 

2.2 

2.4(4) 

0.054 

0.00 

0.63 

0.59 

0.19 

N 

VLASS-D 

2.3 

1.3(4) 

0.887 

0.00 

0.75 

0.75 

0.00 

N 

0.30 

2.8(4) 

0.256 

0.00 

0.57 

0.55 

0.23 

N 

GMB 

0.018 

236(1096) 

0.018 

0.00 

0.88 

0.87 

0.89 

N 

0.0 

727(1096) 

0.004 

0.00 

0.63 

0.63 

0.74 

N 

GMB-Sl 

0.10 

24(110) 

0.033 

0.00 

0.89 

0.87 

0.88 

N 

0.06 

72ni0) 

0.007 

0.00 

0.62 

0.62 

0.63 

N 

CMB-S2 

0.44 

3(11) 

0.054 

0.00 

0.87 

0.86 

0.16 

N 

0.03 

8 (11) 

0.013 

0.00 

0.61 

0.60 

0.58 

N 


Note. — “ Total number of detected events. ^ Average number of detected epochs (total epochs). Average redshift of detected events. Fraction 
of total detected events discarded for lack of variability. Fraction of events with measured rise times. ^ Fraction of events with measured fall times. ® 
Fraction of events with a measured peak. ^ Is galaxy contamination an issue? (Yes, No, Maybe). 


TABLE 4 
LLGRBs 


Survey 

N 

n 

z 

/v 

/rise 

/fall 

/peak 

gal? 

LOEAR 

0.0 

17(37) 

0.001 

0.00 

0.82 

0.86 

0.81 

Y 

LOEAR-E 

0.0 

17(37) 

0.002 

0.00 

0.81 

0.85 

0.80 

Y 

SKA-Low 

18 

6(13) 

0.036 

0.00 

0.81 

0.83 

0.68 

Y 

VAST-W 

0.40 

140(1096) 

0.005 

0.00 

0.96 

0.96 

0.96 

Y 

VAST-W-S 

0.36 

6(37) 

0.011 

0.00 

0.95 

0.95 

0.30 

Y 

VAST-D 

0.18 

1.2(4) 

0.012 

0.00 

0.73 

0.76 

0.00 

Y 

VAST-D-SE 

0.0 

141(1096) 

0.016 

0.00 

0.96 

0.97 

0.96 

Y 

VAST-D-SE-S 

0.03 

6(37) 

0.031 

0.00 

0.95 

0.94 

0.29 

Y 

SKA 

6.6 

2(13) 

0.029 

0.00 

0.90 

0.93 

0.02 

M 

SKA-E 

64 

6(37) 

0.040 

0.00 

0.94 

0.95 

0.30 

M 

VLASS-W 

0.060 

1.1(4) 

0.010 

0.00 

0.77 

0.72 

0.00 

M 

VLASS-D 

0.012 

1.1(4) 

0.053 

0.00 

0.74 

0.77 

0.00 

Y 

GMB 

0.04 

4(1096) 

0.005 

0.00 

1.00 

1.00 

0.00 

Y 

GMB-Sl 

0.034 

1 (110) 

0.008 

0.00 

1.00 

1.00 

0.00 

Y 

GMB-S2 

0.016 

1 (11) 

0.012 

0.00 

0.89 

0.89 

0.00 

Y 


Note. 


Columns are as in Table 3. 
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TABLE 5 
SGRBs 


Survey 

N 

n 

z 

/v 

/rise 

/fall 

^peak 

gal? 

N 

n 

z 

/v 

frise 

/fall 

/peak 

gal? 




^obs 

= 0.2 






^obs 

= 0.4 






LOFAR 

0.0 

12(37) 

0.003 

0.00 

0.83 

0.83 

0.82 

N 

0.0 

16(37) 

0.002 

0.00 

0.77 

0.76 

0.76 

M 

LOFAR-E 

0.0 

12(37) 

0.010 

0.00 

0.83 

0.85 

0.83 

N 

0.006 

16(37) 

0.008 

0.00 

0.75 

0.76 

0.73 

M 

SKA-Low 

22 

4(13) 

0.173 

0.00 

0.81 

0.83 

0.53 

N 

28 

6(13) 

0.125 

0.00 

0.73 

0.75 

0.69 

M 

VAST-W 

0.14 

153(1096) 

0.034 

0.00 

0.91 

0.92 

0.92 

N 

0.06 

270(1096) 

0.018 

0.00 

0.85 

0.86 

0.86 

N 

VAST-W-S 

1.8 

5(37) 

0.082 

0.00 

0.91 

0.91 

0.83 

N 

0.90 

9(37) 

0.043 

0.00 

0.86 

0.83 

0.84 

N 

VAST-D 

1.5 

1.0(4) 

0.098 

0.00 

0.74 

0.74) 

0.00 

N 

1.1 

1.1(4) 

0.051 

0.00 

0.75 

0.75 

0.00 

N 

VAST-D-SF 

0.016 

160(1096) 

0.110 

0.00 

0.91 

0.92 

0.92 

N 

0.006 

276(1096) 

0.057 

0.00 

0.85 

0.85 

0.86 

N 

VAST-D-SF-S 

0.20 

6(37) 

0.255 

0.00 

0.91 

0.90 

0.89 

N 

0.078 

10(37) 

0.129 

0.00 

0.85 

0.85 

0.83 

N 

SKA 

66 

2.1(13) 

0.258 

0.00 

0.89 

0.90 

0.04 

N 

30 

3(13) 

0.135 

0.00 

0.84 

0.84 

0.45 

N 

SKA-E 

420 

6.1(37) 

0.335 

0.00 

0.90 

0.91 

0.89 

N 

170 

10(37) 

0.169 

0.00 

0.85 

0.84 

0.83 

N 

VLASS-W 

1.2 

1.0(4) 

0.097 

0.00 

0.77 

0.75 

0.00 

N 

0.64 

1.0(4) 

0.044 

0.00 

0.77 

0.76 

0.00 

N 

VLASS-D 

0.6 

1.0(4) 

0.707 

0.00 

0.77 

0.74 

0.00 

N 

0.18 

1.0( 4) 

0.264 

0.00 

0.76 

0.74 

0.00 

N 

CMB 

0.0 

150(1096) 

0.004 

0.00 

0.92 

0.92 

0.92 

N 

0.0 

270(1096) 

0.002 

0.00 

0.86 

0.87 

0.87 

N 

CMB-Sl 

0.0 

15(110) 

0.007 

0.00 

0.92 

0.92 

0.92 

N 

0.0 

27(110) 

0.004 

0.00 

0.86 

0.85 

0.86 

N 

CMB-S2 

0.006 

2(11) 

0.012 

0.00 

0.89 

0.91 

0.00 

N 

0.0 

3(11) 

0.007 

0.00 

0.84 

0.84 

0.29 

N 




^obs 

= 0.8 






^obs 

= 1.57 






LObAR 

0.0 

25(37) 

0.001 

0.01 

0.58 

0.61 

0.62 

M 

0.0 

33(37) 

0.001 

0.39 

0.57 

0.43 

0.19 

M 

LOFAR-E 

0.006 

25(37) 

0.004 

0.01 

0.60 

0.59 

0.60 

M 

0.0 

33(37) 

0.002 

0.38 

0.50 

0.50 

0.15 

M 

SKA-Low 

28 

9(13) 

0.070 

0.01 

0.57 

0.58 

0.59 

M 

7.0 

12(13) 

0.029 

0.41 

0.55 

0.45 

0.14 

M 

VAST-W 

0.008 

587(1096) 

0.005 

0.00 

0.69 

0.70 

0.76 

M 

0.0 

960(1096) 

0.001 

0.38 

0.55 

0.45 

0.53 

M 

VAST-W-S 

0.11 

20(37) 

0.012 

0.00 

0.69 

0.70 

0.70 

M 

0.006 

32(37) 

0.003 

0.41 

0.55 

0.45 

0.22 

M 

VAST-D 

0.24 

2.0(4) 

0.016 

0.00 

0.67 

0.68 

0.09 

M 

0.016 

3.4(4) 

0.004 

0.41 

0.53 

0.47 

0.11 

M 

VAST-D-SF 

0.0 

601(1096) 

0.017 

0.00 

0.70 

0.72 

0.76 

M 

0.0 

950(1096) 

0.004 

0.41 

0.53 

0.47 

0.49 

M 

VAST-D-SF-S 

0.010 

20(37) 

0.036 

0.00 

0.70 

0.66 

0.66 

M 

0.0 

32(37) 

0.008 

0.43 

0.56 

0.44 

0.23 

M 

SKA 

3.5 

7(13) 

0.038 

0.00 

0.68 

0.69 

0.63 

M 

0.20 

11(13) 

0.009 

0.42 

0.61 

0.39 

0.20 

M 

SKA-E 

20 

20(37) 

0.048 

0.00 

0.69 

0.70 

0.69 

M 

1.0 

33(37) 

0.011 

0.44 

0.55 

0.45 

0.26 

M 

VLASS-W 

0.072 

1.9(4) 

0.011 

0.00 

0.68 

0.71 

0.04 

M 

0.002 

3.4(4) 

0.002 

0.41 

0.59 

0.41 

0.15 

M 

VLASS-D 

0.014 

1.9(4) 

0.062 

0.00 

0.71 

0.72 

0.06 

M 

0.0 

3.4(4) 

0.012 

0.41 

0.57 

0.43 

0.15 

M 

CMB 

0.0 

591(1096) 

0.001 

0.00 

0.72 

0.71 

0.77 

N 

0.0 

825(1096) 

0.000 

0.00 

0.64 

0.80 

0.78 

N 

CMB-Sl 

0.0 

60(110) 

0.001 

0.00 

0.68 

0.73 

0.72 

N 

0.0 

82(110) 

0.001 

0.00 

0.65 

0.78 

0.75 

N 

CMB-S2 

0.0 

6(11) 

0.002 

0.00 

0.67 

0.67 

0.61 

N 

0.0 

8(11) 

0.001 

0.06 

0.65 

0.64 

0.65 

N 


Note. — Columns are as in Table 3. 


TABLE 6 

NS-NS Mergers with Prompt Black Hole Formation 


Survey 

N 

n 

Z 

/v 

/rise 

/fall 

/peak 

gal? 

LOFAR 

0.0 

36(37) 

0.001 

0.82 

1.00 

0.00 

0.23 

M 

LOFAR-E 

0.016 

37(37) 

0.003 

0.81 

1.00 

0.00 

0.21 

M 

SKA-Low 

48 

13^3) 

0.046 

0.84 

1.00 

0.00 

0.04 

M 

VAST-W 

0.0 

1095(1096) 

0.002 

0.87 

1.00 

0.00 

0.33 

M 

VAST-W-S 

0.060 

37(37) 

0.005 

0.86 

1.00 

0.00 

0.15 

M 

VAST-D 

0.13 

3.9(4) 

0.008 

0.86 

1.00 

0.00 

0.00 

M 

VAST-D-SF 

0.0 

1095(1096) 

0.007 

0.80 

1.00 

0.00 

0.30 

M 

VAST-D-SF-S 

0.006 

37(37) 

0.015 

0.84 

1.00 

0.00 

0.13 

M 

SKA 

1.9 

13^3) 

0.018 

0.84 

1.00 

0.00 

0.08 

M 

SKA-E 

12 

37(37) 

0.020 

0.82 

1.00 

0.00 

0.09 

M 

VLASS-W 

0.026 

4.0(4) 

0.004 

0.83 

1.00 

0.00 

0.00 

N 

VLASS-D 

0.0 

3.9(4) 

0.021 

0.86 

1.00 

0.00 

0.00 

N 

CMB 

0.0 

1096(1096) 

0.000 

0.80 

1.00 

0.00 

0.48 

N 

CMB-Sl 

0.0 

iionio) 

0.001 

0.83 

1.00 

0.00 

0.20 

N 

CMB-S2 

0.0 

11(11) 

0.001 

0.86 

1.00 

0.00 

0.02 

N 


Note. 


Columns are as in Table 3. 
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TABLE 7 

NS-NS Mergers with a Stable Magnetar Remnant 


Survey 

N 

n 

Z 

/v 

/rise 

/fall 

/peak 

gal? 

LOFAR 

0.0 

33(37) 

0.038 

0.76 

1.00 

0.00 

0.00 

N 

LOFAR-E 

0.90 

33(37) 

0.131 

0.76 

1.00 

0.00 

0.14 

N 

SKA-Low 

1160 

12(13) 

1.425 

0.86 

1.00 

0.00 

0.15 

N 

VAST-W 

1.2 

870(1096) 

0.140 

0.58 

0.95 

0.05 

0.73 

N 

VAST-W-S 

12 

31(37) 

0.296 

0.65 

1.00 

0.00 

0.61 

N 

VAST-D 

22 

3.2(4) 

0.371 

0.70 

1.00 

0.00 

0.31 

N 

VAST-D-SF 

0.080 

945(1096) 

0.388 

0.65 

0.99 

0.01 

0.82 

N 

VAST-D-SF-S 

0.72 

33(37) 

0.755 

0.72 

1.00 

0.00 

0.53 

N 

SKA 

220 

12(13) 

0.805 

0.76 

1.00 

0.00 

0.50 

N 

SKA-E 

1020 

34(37) 

0.887 

0.75 

1.00 

0.00 

0.40 

N 

VLASS-W 

5.0 

3.1(4) 

0.227 

0.63 

0.98 

0.02 

0.39 

N 

VLASS-D 

0.38 

3.5(4) 

0.948 

0.76 

1.00 

0.00 

0.16 

N 

CMB 

0.006 

917(1096) 

0.024 

0.47 

0.77 

0.23 

0.72 

N 

CMB-Sl 

0.032 

92(110) 

0.042 

0.47 

0.79 

0.21 

0.67 

N 

CMB-S2 

0.13 

9(11) 

0.070 

0.60 

1.00 

0.00 

0.62 

N 

Note. — Columns are 

as in Table 3. 

TABLE 8 
TDEs 






Survey 

N 

n 

z 

/v 

/rise 

/fall 

/peak 

gal? 

N 

n 

z 

/v 

/rise 

/fall 

/peak 

gal? 



SwJ1644+57 (On 

Axis) 





Off-Axis 






LOFAR 

0.0 

27(37) 

0.013 

0.00 

0.62 

0.64 

0.62 

N 

0.0 

29(37) 

0.021 

0.90 

1.00 

0.00 

0.28 

N 

LOFAR-E 

0.0 

27(37) 

0.042 

0.00 

0.59 

0.63 

0.59 

N 

0.20 

30(37) 

0.072 

0.91 

1.00 

0.00 

0.28 

N 

SKA-Low 

1.9 

11(13) 

0.744 

0.14 

0.48 

0.52 

0.35 

N 

280 

11(13) 

1.271 

0.84 

1.00 

0.00 

0.20 

N 

VAST-W 

0.062 

871(1096) 

0.175 

0.02 

0.50 

0.51 

0.52 

N 

2.0 

902(1096) 

0.129 

0.28 

0.64 

0.42 

0.63 

N 

VAST-W-S 

0.52 

29(37) 

0.424 

0.27 

0.58 

0.42 

0.45 

N 

13 

31(37) 

0.275 

0.39 

0.66 

0.35 

0.63 

N 

VAST-D 

0.82 

2.9(4) 

0.561 

0.46 

0.71 

0.29 

0.31 

N 

20 

3.2(4) 

0.332 

0.54 

0.74 

0.26 

0.38 

N 

VAST-D-SF 

0.0 

289(366) 

0.595 

0.82 

1.00 

0.00 

0.07 

N 

0.018 

359(366) 

0.320 

0.82 

1.00 

0.00 

0.23 

N 

VAST-D-SF-S 

0.0 

10(13) 

1.542 

0.90 

1.00 

0.00 

0.00 

N 

0.070 

13(13) 

0.507 

0.88 

1.00 

0.00 

0.08 

N 

SKA 

2.6 

8(13) 

1.528 

0.76 

0.95 

0.05 

0.41 

N 

72 

11(13) 

0.624 

0.72 

0.90 

0.10 

0.50 

N 

SKA-E 

8.4 

22(37) 

1.641 

0.80 

0.94 

0.06 

0.40 

N 

240 

31 (37) 

0.713 

0.71 

0.90 

0.10 

0.53 

N 

VLASS-W 

0.80 

2.3(4) 

0.719 

0.54 

0.94 

0.06 

0.25 

N 

6.0 

3.2(4) 

0.202 

0.38 

0.68 

0.33 

0.34 

N 

VLASS-D 

0.0 

2.1(4) 

1.697 

0.74 

0.98 

0.02 

0.27 

N 

0.11 

3.3(4) 

0.742 

0.74 

0.91 

0.09 

0.32 

N 

CMB 

0.22 

580(1096) 

0.359 

0.00 

0.82 

0.92 

0.82 

N 

0.012 

903(1096) 

0.021 

0.18 

0.65 

0.52) 

0.66 

N 

CMB-Sl 

0.86 

60(110) 

0.649 

0.00 

0.80 

0.93 

0.79 

N 

0.060 

89(110) 

0.038 

0.22 

0.64 

0.47 

0.64 

N 

CMB-S2 

2.0 

7(11) 

1.036 

0.00 

0.69 

0.87 

0.02 

N 

0.24 

9(11) 

0.065 

0.28 

0.59 

0.42 

0.57 

N 


Note. — Columns are as in Table 3. 


TABLE 9 

Type Ib/c SNe 





Survey 

N 

n 

z 

/v 

/rise 

/fall 

/peak 

gal? 

LOFAR 

0.0 

30(37) 

0.001 

0.07 

0.59 

0.54 

0.60 

Y 

LOFAR-E 

0.0 

30(37) 

0.002 

0.12 

0.60 

0.54 

0.62 

Y 

SKA-Low 

1.6 

11(13) 

0.007 

0.13 

0.60 

0.53 

0.59 

Y 

VAST-W 

0.0 

386(1096) 

0.001 

0.00 

0.88 

0.91 

0.88 

Y 

VAST-W-S 

0.032 

13(37) 

0.002 

0.00 

0.88 

0.90 

0.86 

Y 

VAST-D 

0.036 

1.8(4) 

0.002 

0.00 

0.71 

0.74 

0.01 

Y 

VAST-D-SF 

0.0 

390(1096) 

0.003 

0.00 

0.88 

0.91 

0.88 

Y 

VAST-D-SF-S 

0.0 

14(37) 

0.006 

0.00 

0.87 

0.92 

0.87 

Y 

SKA 

0.85 

5( 13) 

0.006 

0.00 

0.84 

0.88 

0.27 

M 

SKA-E 

1.6 

13( 37) 

0.008 

0.00 

0.87 

0.90 

0.85 

M 

VLASS-W 

0.010 

1.4(4) 

0.002 

0.00 

0.71 

0.78 

0.00 

M 

VLASS-D 

0.0 

1.5(4) 

0.010 

0.00 

0.72 

0.78 

0.00 

Y 

CMB 

0.008 

31(1096) 

0.001 

0.00 

1.00 

1.00 

0.22 

Y 

CMB-Sl 

0.008 

12(110) 

0.002 

0.00 

0.97 

0.98 

0.19 

Y 

CMB-S2 

0.006 

3(11) 

0.003 

0.06 

0.65 

0.93 

0.00 

Y 


Note. 


Columns are as in Table 3. 
















Extragalactic Transients in the Era of Wide-Field Radio Surveys. I 


23 


APPENDIX 


LIMITS FROM PREVIOUS SURVEYS 

Previous surveys are placed in Figure 3 according to their source detection thresholds'^, and their constraints 

of Naii_sky Below we give sensitivity numbers for these surveys. In the vast majority of cases no sources are detected, 
and we take the 95 per cent Poisson credible interval on Naii-sky given this fact. That limit is 3 sources in the survey 
area (Frail et al. 2012), resulting in a lower limit of 


.^all-sky < 123800 



(Al) 


where fig is the effective survey area, which is a function of the transient evolutionary time scale tdur because repeated 
visits to a field do not add more information if the time between them is much shorter than tdur (Macquart 2014). We 
make the approximation that when new epochs are obtained at a cadence Tggp ^ tdun the effective survey area is 

Og = (A2) 

^dur 


where llfp is footprint area of the survey and Tbi is the total time baseline of the survey; the variability loss factor 
V « tdur/7bi- A survey with a duration small compared to tdur has fig —> 0, since separate observations result in zero 
independent measurements in this case. 

Based on the analysis in section 2, here we assume a 100-day timescale in all cases. This means that in many 
cases our adopted values of fig are significantly smaller than the limiting case of the “two-epoch effective survey area” 
(Bower et al. 2007): 

= f^fp(iVep - 1), (A3) 

where Ngp is the number of survey epochs, which applies only when Tgep ^ tdur- In the case of ASKAP-Wide, fig as 
estimated this way is more than two orders of magnitude too large. We note that because both Equations A2 and A3 
have fig oc flfp oc t“^, however, wide and shallow surveys are preferred for Euclidean sources regardless of the relevant 
time scales. 


VLSS (150 MHz) 

Jaeger et al. (2012a) summarize ongoing work to perform a systematic search for transients in the VLA Low Frequency 
Sky Survey (VLSS; Lane et al. 2012). Jaeger et al. (2012b) report one detection, which occurred in one epoch of a 
6-epoch search for transients in a 6.5 deg^ field of view. The detection at 5.5(T (2.1 mJy) occurred in the middle of 
the campaign, which lasted only 3 months. The source appears to rise in flux by a factor of two on a 6-hr time scale 
in the detection epoch. Upper limits for the other epochs are not provided by Jaeger et al. (2012b). 

As the time scale of this survey is about the same as our 100-day choice, we set the detection density to less than 
one source per 6.5 deg^, corresponding to IVaii-sky < 6400. Note that this is well above the predictions for most of our 
transient classes. The rapid variability observed in the discovery epoch does not necessarily indicate unrealistically fast 
timescales, as this may be symptomatic of the low significance of the detection. To summarize, we take = 2.1 mJy 
and iVaii-sky < 6400. 


FIRST (1.4 GHz) 

No transients are reported in comparison of the zero-declination strip of the FIRST survey by de Vries et al. (2004), 
which was visited in two epochs separated by 7 years. The quoted area is 120.2 deg^ (their section 2) and the source 
flux density limit is 2 mJy (their section 2.1.1). We take = 2 mJy and 7Vaii_sky < 1030. 

ATATS/NVSS (1-4 GHz) 

Croft et al. (2010b) compare the Allen Telescope Array Twenty-centimeter Survey (ATATS) with NVSS, reporting 
no transients with a source flux density limit of 40 mJy. As the time baseline between the two surveys is ^ 100 days, 
the effective area is the ATATS footprint of 690 deg^. We thus take = 40 mJy and A^aii-sky < 179. 

FIRST/NVSS (1.4 GHz) 

Levinson et al. (2002) searched for transients in the overlap of the FIRST and NVSS surveys. The effective search 
area is 2400 deg^, and the characteristic limit is 6 mJy. Since FIRST and NVSS were conducted quasi-simultaneously, 
that the time baseline between the two epochs varied between ^0.5 and 4.5 yr. Gal-Yam et al. (2006) followed up the 
transient candidates from Levinson et al. (2002), concluding that two were real: one radio supernovae and a second 
object without an apparent host galaxy. Based on our calculations, the RSN detection was a highly improbable result. 
To summarize, we take Fb.iim = 6 mJy and A"aU-sky = 34. 

Note that in these works it is not necessarily true that = lOcr, as we have adopted. 
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VLA Low-Galactic-Latitude Survey (5 GHz) 

Ofek et al. (2011) report a survey with a footprint of 2.66 deg^ comprising 16 epochs. The spacing of the epochs 
(their Table 3) is such that only 2 of them are independent on 100-day time scales, so we adopt A = 5.32 deg^, 
significantly smaller than the naive effective area of 41 deg^. An Scr source was reported that was present in the very 
first epoch and became undetectable 3 days later. Since this evolution appears unrealistically rapid, this source was 
probably spurious even though it passed various quality checks, a conclusion shared by other groups (Frail et al. 2012). 
Ofek et al. (2011) report a representative flux density limit of 1.8 mJy (their section 9.1), so we take = 1.8 mJy 

and TVaii-sky = 23270. We note that although we present the limit from this survey with our 3 GHz results, the survey 
was conducted at 5 GHz. 


VLA Galibrator Field Survey (5 & 8 GHz) 

Bower et al. (2007) report the results of a 944-epoch survey for transient sources with archival data from the VLA 
spanning 22 years with a typical epoch separation of 7 days. We update the results to account for the improved 
re-analysis of Frail et al. (2012). At a 100 day cadence, the 22 years of observations represent only 80 independent 
epochs, so the effective area is reduced by a factor of ^ 10 compared to a naive estimate. Observations were performed 
at 4.9 and 8.4 GHz, and again we show this survey with our 3 GHz results. 

Frail et al. (2012) conclude that there is only one transient in the data set, with a flux density of 0.6 mJy at 4.9 GHz 
and a significance of 5.8cr. It only appears in one epoch. As with the Jaeger et al. (2012b) transient, the detected event 
is too fast and implies a rate far above our predictions, but we treat it as real. Bower et al. (2007) report effective 
areas as a function of the flux density limit and the characteristic time scale (their Tables 5-7). We employ the limit 
closest to the detected transient flux density of 0.56 mJy. On a 2-month time scale, their Table 6 gives an effective 
area of 4.97 deg^ with 120 epochs. Scaling by 80/120 (i.e., converting Vlg to use the definition of Equation A2 rather 
than Equation A3) gives a corrected effective area of 3.3 deg^. We thus take Fj, um = 0.56 mJy and IVaii-sky = 12500. 
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